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Giant radio relics in galaxy clusters: reacceleration of fossil 
relativistic electrons? 
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ABSTRACT 

Many bright radio relics in the outskirts of galaxy clusters have low inferred Mach 
numbers, defying expectations from shock acceleration theory and heliospheric ob- 
servations that the injection efficiency of relativistic particles plummets at low Mach 
numbers. With a suite of cosmological simulations, we follow the diffusive shock ac- 
celeration as well as radiative and Coulomb cooling of cosmic ray electrons during 
the assembly of a cluster. We find a substantial population of fossil electrons. When 
reaccelerated at a shock (through diffusive shock acceleration), they are competitive 
with direct injection at strong shocks and overwhelmingly dominate by many orders of 
magnitude at weak shocks, 3, which are the vast majority at the cluster periph- 

ery. Their relative importance depends on cooling physics and is robust to the shock 
acceleration model used. While the abundance of fossils can vary by a factor of ~ 10, 
the typical reaccelerated fossil population has radio brightness in excellent agreement 
with observations. Fossil electrons with 1 r$ 7 ^ 100 (10 ^ ^ ^ 10^) provide the main 
seeds for reacceleration at strong (weak) shocks; we show that these are well-resolved 
by our simulation. We construct a simple self-similar analytic model which assumes 
steady recent injection and cooling. It agrees well with our simulations, allowing rapid 
estimates and physical insight into the shape of the distribution function. We predict 
that LOFAR should find many more bright steep-spectrum radio relics, which are 
inconsistent with direct injection. A failure to take fossil cosmic ray electrons into 
account will lead to erroneous conclusions about the nature of particle acceleration 
at weak shocks; they arise from well-understood physical processes and cannot be 
ignored. 

Key words: magnetic fields, cosmic rays, radiation mechanisms: non-thermal, ele- 
mentary particles, galaxies: cluster: general 



1 INTRODUCTION 

Diffuse radio emission in clusters falls into two broad 
clas ses: smooth, centra lly located and unpolarized radio ha- 
los (|Ferrari et al.|[200^ . and references therein), and elon- 
gated, significantly polarized and stee p spectrum radio relics 
which are seen at cluster outskirts jKempner et al.l |2004| . 
and references therein). The radio halos come in two dis- 
tinct classes: mini-halos in cool core clusters and giant ha- 
los associated with merging clusters. Similar to the ra- 
dio halos, the radio relics are thought to come in dis- 
tinct classes which are all associated with merging clus- 
ters: fossil radio lobes blown in the past by active galac- 
tic nuclei (AGN), which may have been re-energized by 
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compression (radio phoenix: lEnfilin fc Gopal-Krishnall200ll ). 
and those produced by direct diffusive par ticle accelera- 
tion at accretion/merger shocks (radio gischt: lEnsslin et al.l 
Il998l : iMiniati et al.ll200il ). Many radio r elics have now been 
seen 50: for r ecent compilations, see lFeretti et al.ll2012l : 
iNuza et al.l I2OI1I ) the number could increase dramatically 
with upcoming low frequency radio surveys proposed for the 
Low Frequency Array (LOFAR), the Westerbork Synthesis 
Radio Telescope (WSRT), and farther in the future. Square 
Kilometer Array (SKA). Using numerical simulations com- 
bined with a semi-analytic model for radio emission cali- 
brated to existing number counts, iNuza et al.l l|201ll ') esti- 
mate that LOFAR and WSRT could discover ~ 2500 relics 
and ~ 900 relics respectively. The time is therefore ripe to 
understand how we can best mine these future surveys. 

In this paper, we focus on radio gischt, which trace 
structure formation shock waves. They therefore can illumi- 
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nate the nature of cosmic accretion/mergers , as well as shock 
ampl if ication of large-scale magnetic fields (Pfrommer et al 



20081 : IPfrommeJ l2008l: iHoeft et all i20Q8.; Batta gha et al. 



200g : ISkillman et al.H201ll '). Perhaps even more importantly, 



they allow us to probe in detail the efficiency of shock accel- 
eration in a diffuse, low Mach number {M ~ 2 — 4) regime, 
far different from the high Mach number regimes probed in 
our Galaxy and supernova remnants. Whilst this remains 
poorly understood, at face value the observations seem to 
suggest an electron acceleration efficiency significantly in ex- 
cess of naive theoretical expectations. 

A recent spectacular example of radio gischt is CIZA 
J2242.8-|-530f, the 'sausage relic' (|van Weeren et al.ll201(3 ). 
a large (~ 2 Mpc long; located ~ f .5 Mpc from the clus- 
ter center) double radio relic system. The post-shock radio 
spectral index was used to infer the particle spectral slope 
and hence the shock compression ratio and Mach number 
{M ~ 4.6), while the decrease in the spectral index (from 
~ 0.6 to ~ 2.0 across the relic's narrow ~ 55 kpc width) to- 
ward the cluster center — spectral aging due to synchrotron 
and inverse Compton losses — was used to infer magnetic 
field strengths of ~ 5 ^G. The strong (~ 50 — 60 per cent) 
polarization can be attributed to magnetic field frozen into 
the compressed intracluster medium (ICM), which has been 
aligned parallel to the shock. The properties of this and 
similar systems are distinct from fossil radio plasma, which 
are smaller, have curved, steeper spect ra (due to aging), 
and lobe- or torus-like morpho logy (e.g.. lEnfilin fc Briiggenl 
l2002l : IPfrommer fc Jon3l2011^ ■ The power-law spectral in- 
dex, spectral gradient, and enormous extent clearly support 
a diffusive shock acceleration (DSA) origin. 

However, this then presents a puzzle. Cosmological 
simulations show that while gas initially undergoes strong 
shocks (up to M fO^) accreting onto non- linear struc- 
tures and filaments, shocks in the ICM and cluster out- 
skirts are relatively modest {Mi ~ 1 — 4 ), since the gas 
has already reached su b -keV temperatu r es (IRvu et al.ll2003l: 
Pfromm er et al.l l2006l : ISkillman et"al] l2008l : IVazza et all 



2009, 2oTii r While DSA is efficient in accelerating parti- 
cles in the t hermal Maxw ellian tail at high Mach numbers 
l|Bellll f 978bl : lDrurvlll983bh . and as confirmed in observations 
of supernova remnants (jParizot et al ] |2006l : lRevnolds|[200i ) , 
at lower Mach numbers the efficiency of DSA is known to 
plummet exponentially. Indeed, in the test-particle regime 
where suprathermal particles undergo acceleration via a 
thermal leakage process (which compare well against kinetic 
DSA simulations) the acceleration efficiency for weak shocks 
< 3 is extremely small; the fraction of protons acceler- 
ated is ~ 10~* — f0~^, and cosmic ray proto n (CRp) pres- 
sure is ^ I per cent of the shock ram pressure (|Kang fc Rvul 
[20l3). The acceleration efficiency of electrons at low Mach 
numbers is likely to be far smaller still. The injection prob- 
lem for thermal electrons is already known to be severe at 
high Mach numbers, due to the smaller gyroradius of ther- 
mal electrons: the relative acceleration efficiency of cosmic 
ray e lectrons (CRes) is lower by ~ f0~^ as in the Galaxy 
IjSchlickeiser 2002) or eve n ~ 10"* as in supernova rem- 
nants ( Morlino et al.ll2009l ). These relative efficiencies likely 
plummets further at low Mach numbers, as is also suggested 
by heliospheric observations (see 37|. These considerations 
appear to contradict the appearance of bright radio relics, 
and perhaps suggest that our understanding of DSA at 



low Mach numbers is incom plete (for recent progress, see 
iGargate fc Spitkovsk-\^|2012| ). 

A possible solution is if there is a pre-existing popu- 
lation of CRes with gyroradii comparable or larger than 
that of the shock thickness. In this case, i njection from 
the th ermal pool is no longe r an issue (lEnsslin et ali 



199 81: iMarkevitch et all l2005l: iGiacintucci et al.l 12001 



Kang fc Rvul 1201 ll : iKang et alj Wlj ). DSA has a much 



larger effect on radio emission than adiabatic compression; 
including the downstream magnetic field amplification, syn- 
chrotron emission from fossil ele ctrons could be boo sted by 
a factor ~ 100-1000 for Al ~ 3 l|Kang fc Rvull201ll '). Thus, 
the luminosity function of radio gischt will be strongly mod- 
ified by the presence or absence of a seed relativistic electron 
population, whose existence has never been directly demon- 
strated. Note that secondary CRe which arise from hadronic 
interactions of CRp are not thought to be signifi cant at 
these low densities, though for a dissenting view, see lKeshe3 
(|2010l ). In this paper, we use our exi sting high-resolution 
SPH simulations of CRps in clusters (|Pinzke fc Pfrommer] 
I2O IO ) to infer whether structure formation shocks could gen- 
erate thertQ. Note that DSA operates identically on rela- 
tivistic particles of the same rigidity {R — Pc/Ze), so the 
injected proton and electron spectrum are the same, modulo 
their relative acceleration efficiency, which we calibrate off 
Galactic observations. 

The main difference is that unlike protons, relativistic 
electrons can still undergo significant Coulomb and inverse 
Compton losses in the cluster outskirts. In Fig. [T] we show 
the cooling time of CRes for difference densities and cos- 
mological epochs. CRes of energy ~ 10 GeV that are re- 
sponsible for ~ GHz emission in a fiG field have energy loss 
timescales of ~ lO^yr at z — (and cool even more quickly 
at higher redshift). However, lower energy CRe could po- 
tentially survive to be reaccelerated; note that CRes with 
momentum P = (1 — 100) rric c have tcooi > in in gas with 
rie = 10~^ cm""^ (typical of cluster outskirts) at z = 0. We 
evolve the time-dependent cosmic ray (CR) energy equation 
to track the evolving distribution function of CRe electrons 
under the combined influence of injection and cooling. We 
find that the fossil electron population is substantial, and 
that reacceleration of this population is competitive with 
direct injection for strong shocks, and will overwhelming 
dominate radio emission for weak shocks with Mach num- 
bers Al < 4. Since the latter constitute the vast majority 
of shocks in the cluster periphery, reacceleration is a crucial 
mechanism which cannot be ignored. 

The outline of this paper is as follows. In we de- 
scribe our CR formalism and computational method. In ij3l 
we present results for the fossil electron spectrum, and an 
analytic model for understanding its essential features. In 
[J4l we study how the fossil population is transformed by a 
shock, time-resolution effects, and our principal result, that 
reacceleration dominates direct injection. In iJS] we study 



^ They could also have a non-gravitational origin, such as accel- 
eration by AGNs or SN-driven winds, though the filling factor 
at these large radii is likely to be small and the effect of adia- 
batic cooling during the expansion from the compact interstellar 
medium into the ICM makes their energy fraction likely negligible 
in comparison to those CRes accelerated at structure formation 
shocks. 
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Figure 1. Cooling time of CR electrons for typical conditions in 
the ICM. CRes can be long-lived if they are injected at low density 
and at late times. Electrons with kinetic energy E ~ prriaC^ ~ 
10 MeV should be the most long-lived in cluster outskirts (jic ~ 
10~^cm~3) because this energy range maximizes the Coulomb 
and inverse Compton (IC) cooling times. 

how the fossil spectrum varies with the shock acceleration 
model and between clusters. In [JG] we discuss observational 
consequences, including the brightness of relics as a func- 
tion of Mach number and the relic luminosity function. In 
fJT] we discuss insights gained from heliospheric observations 
on our adopted shock physics as well as possible limitations 
of our numerical method. We conclude in SJS] Three appen- 
dices also spell out various technical points about the cooling 
and injection process. 



2 FOSSIL ELECTRONS: METHOD 

Our simulations adopt a ACDM cosmology with parame- 
ters: = !^DM + Ob = 0.3, Ob = 0.039, Oa = 0.7, h = 
0.7, ris = 1, and as = 0.9. The simulations were carried out 
with an updated and extended version of t he distributed- 
memory parallel Tre eSPH code GADGET-2 (|Springel2005l : 
ISpringel et alll200ll ). Gravitational forces are computed us- 
ing a combination of particle-mesh and tree algorithms. 
Hydrodynamic forces are computed with a variant of the 
smoothed particle hydrodynamics (SPH) algorithm that 
conserves energy and entropy where appropriate, i.e. outside 
of shocked regions (Springcl & Hcrnquist 2002). Our simula- 
tions follow the radiative cooling of the gas, star formation, 
supernova feedback, and a photo-ionizing b ackground (de- 
tails can be found in iPfrommer et al.ll2007l). We model the 
CR physics in a s e lf-consistent w ay "("Pfromme r et al.ll2006l : 
lEnfilin et ai]|2007l : lJubelgafi et al. 2008) and attach a CRp 
distribution function to each SPH fluid element. We include 
adiabatic CRp transport process such as compression and 
rarefaction, and a number of physical source and sink terms 
which modify the CRp pressure of each CRp population sep- 
arately. The main source of injection is diffusive shock accel- 
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Table 1. Cluster sample. 
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M 
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1.7 


1.049 
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1.6 


1.049 
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M 


9.1 X 10^3 


0.93 


1.3 


1.127 
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M 


8.5 X 10^3 


0.91 


1.3 


1.127 


gSd 


PreM 


7.8 X 10^3 


0.88 


1.2 
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Notes: 



(1) The dynamical state has been classified through a combined 
criterion invoking a merger tree study and visual inspection of the 
X-ray brightness maps. The labels for the clusters are: PreM— pre- 
merger (sub-cluster already within the virial radius), M-merger, 
PostM-post-merger (slightly elongated X-ray contours, weak cool 
core (CC) region developing), CC-cool core cluster with extended 
cooling region (smooth X-ray profile). (2) The virial mass and ra- 
dius are related by Ma(^) = ^tt A pcrit(2:)fi^, where A = 200 
denotes a multiple of the critical overdensity Pcrit . (3) The virial 
temperature is kT200 = 0^/200 A'"ip/(2-R20o), where /i is the 
mean molecular weight. (4) Time difi'erence between output snap- 
shots; in units of Myr for g72, and remaining clusters show the 
ratio of the cosmological scale factor, a, between two snapshots 
(roughly corresponding to a time interval At ~ IGyr). 

eration (DSA) at cosmological structure formation shocks, 
while the primary CRp sinks are thermalization by Coulomb 
interactions, and catastrophic losses by hadronization. We 
do not consider CRps injected by supernovae remnants or 
AGN feedback; these sources should be relatively subdomi- 
nant in the cluster outskirts. 

In this paper, we post-process previous cosmologi- 
cal s imulations of 14 galaxy clusters (|Pinzke fc Pfrommeij 
ioiO), adapting them to study radio relics. The properties 
of the cluster sample is listed in Table [T] In our previous 
work, CRps were accelerated through DSA on the fly in the 
simulations, and the CRp distribution was written out into 
snapshots. We now use this information to identify shocks 
and inject CRes. In particular, we; 

• identify all SPH particles that have undergone a shock 
by comparing the CRp distribution function between snap- 
shots 

• inject CRes according to acceleration scheme (described 
in i^TT^ 

• evolve each injected CRe population to a later time tf 
while accounting for losses (Coulomb, inverse Compton, and 
adiabatic) 

• add up the CRe distribution for all SPH particles at 
time tf. this constitutes our fossil spectrum. 

• reaccelerate CRe distribution function at time tf using 
typical values for a merging shock in cluster outskirts 

• calculate the radio synchrotron emission from radio 
relics 
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In this section, we describe our scheme for deriving the 
fossil CRe distribution function at time tf. Reacceleration 
and radio emission are described in ^and |6] respectively. 

2.1 Basic CR Formalism 

The CR electrons and protons are each represented by an 
isotropic one-dimensional power-law distribution function 
given by 



fip) 



dpdV 



(1) 



where p = P/m c, where P is the momentum, m the mass of 
the particle, and c the speed of light. Note that for Lorentz 
factors 7 2> 1, p ~ 7. Additionally, pcut is the momentum 
cutoff, a the spectral index of the CR power-law distribu- 
tion, C the normalization of the distribution function, and 
H{x) is the Heaviside step function. The differential CR 
spectrum can vary spatially and temporally, but for brevity 
we suppress this in our notation. 

The CR number density is given by 



ncR 



f 

Jo 



Apf{p) = 



CpI-C 



(2) 



provided a > 1. The kinetic energy density of the CR pop- 
ulation is 



ecR = 



'dp/(p)T(p) = ^l^ X 
a — 1 



a - 2 3 - a X 

+ Pcut 



(3) 



where T{p) = {y/TTr? ~ 1) m c? is the kinetic energy of 
a particle with momentum p and x — 1/(1 -|-Pcut)- ^x{oL,h) 
denotes the incomplete Beta-function, and q > 2 is assumed. 
The average CR kinetic energy Tor = ecr/^-cr is therefore 



TCF 



2 3 



0^/ 



1+Pcut-l 



(4) 



2.2 CR Injection 

2.2.1 Identifying Shocks 

In our cosmological simulations, we model the CR proton 
(CRp) distribution function as a superposition of CRp pop- 
ulations, each determined by equation ([l|, but with a differ- 
ent spectral index {a — {2.1,2.3,2.5,2.7,2.9}), momentum 
cut-off Pcut, and normalization Cp (where the subscript de- 
notes the CR proton population) derived from the simula- 
tions. Both Cp = C{x,t) and pcut = p{x,t) are allowed to 
vary spatially and temporally. 

We identify the SPH particles that have experienced a 
shock by calculating the change in the adiabatic invariant of 
the CR normalization, Co,p = {p/ po)~^°'^^^^^Cp, between a 
snapshot at time t and an earlier time t — At: 

ACo,p(f) = Co,p(t)-Co,p(t- At). (5) 



^ The three-dimensional distribution function is /('^' (p) 
/(p)/(4^rp2). 



Here the time between snapshots is denoted by At. It is 
shown in Table [1] for each simulated cluster. We resolve 
cooling and injection in two of the simulated clusters (g72a 
and g72b) on At ~100 Myr timescales, while the remain- 
ing clusters in our sample are resolved on longer timescales 
At 1 Gyr. If the condition ACo,p(t) > is fulfilled, then 
the CRps have experience at least one shock within the time 
At, since the cooling time of CRps on the cluster outskirts 
is sufficiently long that they are essentially adiabatic. 

2.2.2 Diffusive Shock Acceleration 

Here we introduce the basic framework of DSA which we 
use. In our discussions {po,uq) and {p2,U2) refer to up- 
stream and downstream densities/velocities respectively. In 
the thermal leakage model for DSA (e.g., Jones fc Kand 
ll993l : lBerezhko et al ] | 19941 : iKang fc Jonesll 19951 ). only parti- 
cles in the exponential tail of the Maxwellian thermal distri- 
bution will be able to cross the shock upstream to undergo 
acceleration. The threshold momentum is: 



Pinj — 2;injPth 



2fcT2 



(6) 



where typically a;inj ~ 3.5 — 4. We adopt a fit to Monte Carlo 
simulations of the thermal leakage process (|Kang fc Rvul 
l201i): 



1.17 



U2 
Pth C 



1 + 



1.07 



M 

3 



(7) 



where M is the Mach number of the shock, ea = Bo/B±, 
Bo is the amplitude of the downstream MHD wave turbu- 
lence, and B± is the magnetic field along the shock nor- 
mal. The physical range of ea is quite uncertain due to 
complex plasma interactions, although both plasma hybrid 
simulations and theo ry suggest that 0.25 ;$ ea ;$ 0.35 
(jMalkov fc Volklll998l ). In this paper, we adopt ea = 0.23, 
which corresponds to a conservative maximum energy ac- 
celeration efficiency for electrons (equation I15|l of ~ 1 per 
cent. We show in Appendix[B]how our results vary with this 
parameter. 

The thermal post-shock distribution is a Maxwellian: 

3/2 



.fth (p) = 47t nth 



p exp 



2 2 
mc p 

' 2fcT2 



(8) 



where the gas number density (nth) and temperature {T2) 
are derived from the shock jump conditions. In the test par- 
ticle regime, the CR power-law attaches smoothly onto the 
thermal post-shock distribution, at Xinj, which is the only 
free parameter: 



,fcR,lin(p) = ./■th(pinj) ( — 
\Pinj 



H{p~Pini). (9) 



Fixing the normalization of the injected CR spectrum by 
this continuity condition automatically determines the nor- 
malization constant Cinj . The slope of the injected CR spec- 
trum is: 



ftinj = 



r + 2 



(7ad + l)M^ 



r - 1 (7ad - l)M^ + 2 
where 7ad = 5/3 is the adiabatic index and 
r — ^ — !^ 

pO U2 



(10) 



(11) 
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Figure 2. Acceleration efficiency as a function of shock strength 
for the spectral index renormalization model ('M-aR', see i|2.2.3l 
for details). Those are compared to the electron acceleration effi- 
ciency in various models: 'M-qR' (red dashed), the test particle 
acceleration model ('M-testp'), and a model with constant accel- 
eration efficiency that is widely used in the literature ('M-const'). 
The thick lines represent shocks in a hot medium with a post- 
shock temperature kT2 = 10 keV, and the thin lines represent 
shocks in colder media with kT2 =0.1 keV. 



denote s the shock compression ratio (|Beil Il978al lbl: iDrurvl 
Il983al ). We assume here that the upstream Alfven Mach 
number Ma = U2/va ^ 1, so magnetic fields are dynami- 
cally unimportant. The number density of injected CR par- 
ticles is given by 



AncR.lin = / dp/cR,lin(p) = /th(Pinj) ■ 



Pii 



(12) 

'0 "inj ~ ^ 

This enables us to infer the particle injection efficiency, 
which is the fraction of downstream thermal gas particles 
which experience diffusive shock acceleration, 

„3 



??CR,lin 



AncR.lin 



1 



(13) 



nth '^inj 
The particle injection efficiency is a strong function of a^inj 
that depends on both A4 and eb (for instance, it changes 
by more than an order of magnitude for eb = 0.25 — 0.30 
at = 3. We discuss this further in Appendix [B|. The 
energy density of CRs that are injected and accelerated at 
the shock (neglecting the CR back reaction on the shock) is 
given by 

AecR.i in — ^CR,lin (A^) rcR(A1,Pinj) nth(T2), (14) 
and the CR energy injection and acceleration efficiency is: 
AecR.iiii 



Clin 



where Asdis 



eth2 



£thor^="^. (15) 



ASdiss 

The dissipated energy density in the downstream regime, 
Asdiss, is given by the difference of the thermal energy den- 
sities in the pre- and post-shock regimes, corrected for the 
adiabatic energy increase due to gas compression. 

2.2.3 Models for the CR Spectrum 

The above formalism for diffusive shock acceleration de- 
scribes the standard test particle scenario in which the 



relativistic particles have no influence on the shock struc- 
ture. However, these equations are no longer valid at high 
Mach numbers, when they predict that the fraction of 
energy which goes into relativistic particles (as in equa- 
tion [15]) can reach or exceed 100 per cent. By contrast, nu- 
merical studies of shock acceleration of ions suggest that 
Clin has a up per energy injectio n efficiency limit set by 
Cmax — 0.1 (jKang fc Rvul 120131 ) . This limit arises from 
nonlinear effects, due to the back-reaction of the acceler- 
ated particles upon the sh ock tei chlei 19791: Drurv &; Voelkl 
1981; Axford et al.l Il982l: Malk ov et al.l l2000l : iBlasil |2002| : 
iBlasi et al.ll2005l : lKang fc Jonesll2005l ). 

While non-linear models of particle acceleration cer- 
tainly exist, for our purposes they introduce needless com- 
plication: as we shall see, due to the strong cooling processes 
at play, we are not strongly sensitive to detailed features of 
the injected spectrum. We therefore employ simple modifi- 
cations of the test-particle picture to ensure energy conser- 
vation is not violated. We require that the energy injection 
efficiency Cinj obeys an upper bound ("max: 



Cinj = 1 - exp - 



Clin 



Cn 



Cn 



(16) 



where we set Cmax = 0.1 for the CRp population, and 
Cmax = 0.01 for the CRe population (equivalently, we can 
set Cinj ~ min(Ciin, Cmax)) □• If the injected one-dimensional 
CRe distribution function at each time t is: 



/inj,CRo(p, t) = ACc P "'"'H{p - Pent) 



(17) 



then there are 3 effective time-dependent variables 
ACc, Oinj, Pent to Satisfy the twin constraints of injected 
number density (equation I13|l and energy density (equa- 
tion [14] and I16|l . We consider two simple models. Model 
'M-aR' is motivated by models of non-linear shock accel- 
eration, and varies (ACo,ainj) (keeping pcnt ~ Pinj fixed), 
while Model 'M-testp' varies (ACe,Pcut), keeping the test 
particle slope ainj constant. We also account for the differ- 
ent acceleration efficiencies of protons and electrons (due to 
their different masses, which give rise to different gyroradii 
in the non-relativistic regime) by assuming the CRe distri- 
bution function mimics the CRp distribution function, but 
with a lower normalization. The prescriptions for these two 
models is: 

• Spectral index renormalization ('M-aR') We 

limit the acceleration efficiency by steepening the spectral 
index of the injected population ainj to asub. The slope 
Qinj affects Cinj via the mean energy per particle, equa- 
tion Q. This is motivated by models of non- linear shock ac- 
celeration where a subshock with a lower com pression ratio 
(and hence steeper spectral index) forms (e.g.. [Ellison et ahl 
I2OOOI ). Given the assumed e_g , we find that for strong shocks 
where a < 2.3 the spectral slope is steepened by a maximum 
of ~ 10 per cent in low temperature regimes (fcT ~ 0.1 keV), 
while the steepening is much smaller for high temperature 
regimes (fcT ~ 10 keV) that are more relevant for clusters. 



^ Note that the subscript 'lin' denotes the injected and accel- 
erated population in the test particle approximation, while 'inj' 
refers to the injected and accelerated population after accounting 
for CR back reaction effects on the shock structure that leads to 
saturation of the accelerated CR energy density. 
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Since pinj remains fixed, so does ncR, and we solve for the 
normalization constant ACc from equations ((2} and (|12|) : 

ACp = r?CR,lin (Qsub - 1) P",-""'' • (18) 

We then relate the injected CRes to the CRps by assuming: 



/e(P) = AVp(a)/p(P), 



(19) 



where fc{P) and fp{P) are the CR distribution functions 
at physical momentum P for the electrons and protons, re- 
spectively. The ratio between the electrons and protons can 
be fixed by requiring equal CRe and CRp number densi- 
ties above a fixed injection energy i?inj, which results in 
Kg/p{a) = X^^""-*/^, where X = m^/ me is the proton 
to electron mass ratio (|Schlickeiseill2002l ). For the value of 
a ~ 2.3 (consistent with the injection spectral index at 
galactic supernova remnants as traced by multi-frequency 
observations), this yields K^/p ~ 0.01, which is what is ob- 
served locally. Note that the appropriate choice of normal- 
ization can be highly uncertain — for instance, supernova 
remnant s give values of K^/^ ~ 10~^ which are significantly 
smaller (jMorlino et al.ll20odh . However, it does not signifi- 
cantly affect our results in this paper, which hinge on the 
scaling of acceleration efficiency with Mach number. 

• Modified test particle ('M-testp') This model pre- 
serves the test particle slope Q^inj and instead modifies 
the cutoff pcut of the power-law distribution, which is no 
longer necessarily equal to pinj (jEnfflin et all |2007 '\ This 
was argued to mimic the ef fect of rapid Coulom b cooling in 
the non-relativistic regime jEnfilin et al.l 120071 ') . though the 
physical motivation for this is less clear (our calculations 
explicitly track Coulomb cooling). Nonetheless, we include 
this model to allow comparison with previous published re- 
sults. We then solve the two equations AecRc(Qiiij, gcut) = 
Cinj.cAediss and AncRc(ainj, gcut) = A rycR.iin nth(72) for Co 
and (jcut. 

We regard the 'M-qR' model as our default model, 
since the steepening of the spectral slope is physically well- 
motivated. In practice, we shall see that these two mod- 
els give almost identical results, as the differences between 
them in the momentum regime important for fossil reac- 
celeration is negligible. We also adopt a straw-man model 
where a constant energy acceleration efficiency independent 
of Mach number is assumed: 

• Constant acceleration efficiency ('M-const') Sim- 
ilar to the 'M-testp' model, we solve for both pcut and ACe 
with a fixed to the test particle slope for a combination of 
ee,ncRe; however, the energy injection efficiency Cinj.c ~ 1 
per cent is assumed to be constant and independent of Mach 
number. 

Unlike the previous two models, we do not regard this model 
as physically realistic. It ignores what we know about the 
strong dependence of acceleration efficiency on Mach num- 
ber, and consequently vastly overestimates the number of 
weak shocks visible if only direct injection operates (see !j6]). 
Nonetheless, because it has been widely adopted in the liter- 
ature, we adopt it as a straw man model to compare against 
previous results. 

We show the acceleration efficiency for all 3 models in 
Fig. [2] The 'M-aR' and 'M-testp' models give very similar 
results. They are comparable to the 'M-const' model for high 
M, but have much lower acceleration efficiencies at low M. 



2.3 CR Electron Cooling 

The CRes cool through synchrotron and inverse Compton 
(IC) emissiorlfl and Coulomb interactions on timescales that 
are relatively short compared to the dynamical timescale of 
a cluster. The finite time resolution of our simulations im- 
ply that we inject electrons at discrete time intervals, rather 
than continuously. We therefore incorporate cooling in our 
simulations by considering the evolution of a power-law of 
spectrum electrons instantaneously injected at time ti and 
evolved forward to a later time tf (for further details, see 
ISarazinI (|l999l ). and Appendix |A)| . If there is no further in- 
jection of CRes, then the number of particles is conserved: 

/inj,CRc(p',i/)dp' = / /inj,CRo(p',t»)dp' . (20) 
' P! J Pi 

The distribution function is then given by 



/inj,CRo(Pf,i/(ti)) = /inj,CRo(Pi,ti) 



dpf 



where (21) 



/mj,CRc(Pi,ti 



/inj,CRc(Pf - Apic - Apcoui,ti)-(22) 



Here Apic and Apcoui represents the shift in momentum 
from Pi to pf of a CRe, due to inverse Compton and Coulomb 
cooling, respectively. It is derived by integrating the loss 
function b(p,t), defined by 



where E denotes the particle energy in units of meC? 
AE/&t represents the loss of energy for each particle. 
The Coulomb losses for CRes are described by 



(23) 
and 



bc{p,t) 



where 



3 ot ric c 
2/?2 



In 



1 



ln(2) ( ^ + - 



1 

+ 2^ 



47 



(24) 



P = p/ x/TTp^ , and 



\/r+i 



Here Wpiasma = \/47r e^ne/rrie is the plasma frequency, and 
rio is the number density of free electrons. The details of 
the Coulomb cooling are given in Appendix]^ An impor- 
tant feature to note is that 6 oc p~^ for non-relativistic (NR) 
electrons, while b ~ const for relativistic (R) electrons, im- 
plying tcooi oc p^ (NR) and tcooi oc p (R). We will see how 



(25) 



this shapes the fossil spectrum in 

The inverse Compton losses are given by 

, , , 4 (tt P^ tt 1 /-. \4 

6ic(p, z) = UcMB = oic,o — ( l + z) 

irUcC p p 

where the energy density of the CMB, C/cmb = 
0.26[eVcm"^](l-|-2)'', and the Thomson cross section, (tt = 



87re /3(mo 



Similarly to Coulomb cooling we derive 



the momentum evolution of a particle subject to IC losses 
through 



-&ic.o(l + z)*dt. 



(26) 



* Hereafter, the initials 'IC should be understood to denote 
both processes. For the purpose of calculating the fossil elec- 
tron distribution function, we ignore synchrotron cooling. The 
CMB has an energy density equivalent magnetic field strength of 
BcMB = 3.24(1 -I- z)^^G, and is generally dominant. 



© 2011 RAS, MNRAS 000.[Tlf23l 



Giant radio relics in galaxy clusters: reacceleration of fossil relativistic electrons? 7 



When a time At — {tf — ti) has elapsed, all energetic CRes 
with a momentum p ^ pic have cooled to a lower momentum 
p < pic- We derive pic by integrating equation (|26[) from 
the redshift Zi = z{ti), where the electrons are injected, to a 
later time Zf = z(tj) where they are evaluated. The details 
of the IC cooling are given in Appendix \K\ 

The total electron spectrum is derived from the sum of 
all individually cooled injected spectra (denoted by summa- 
tion index j) Q starting from the time of injection ti until a 
later time t[, 

/CRe(Pf,t/) = ^/inj,CRo(Pf,t/,tj) • (27) 
j 

It is important to note that our injection time steps 
of At = 100 Myr, 1 Gyr are in fact often longer than the 
cooling time of electrons at both low and high energy ex- 
tremes (see Fig. [T]). This low time resolution can have the 
effect of severely overestimating the impact of cooling. Even 
so, it turns out that our simulations generally give results 
accurate to within a factor ~ 2 for the reaccelerated spec- 
trum, compared to calculations where injection and cooling 
are fully resolved. We discuss this in detail in ii4.2l 



3 FOSSIL ELECTRONS: RESULTS AND 
PHYSICAL INTERPRETATION 

We present our primary results for the relic electron dis- 
tribution function in this section. In ^ we analyze a rep- 
resentative CR electron spectrum in the cluster outskirts, 
(1.3-1.5)7?2oo 0, 

that we obtain from our simulations. We 
focus on the underlying physics which gives rise to the spec- 
trum, and show in H3.2l that a simple analytic model can fit 
our simulation results in the momentum regime resolved by 
our simulations. 



3.1 Representative Fossil Spectrum: Simulation 
Results 

Let us first consider the spectrum from a typical cluster, the 
post-merger cluster g72a, generated using our 'M-aR' model 
(from Table [TJ note that 9 of the 14 simulated clusters are in 
the process of merging or had a recent merger) . We focus on 
a particular case because the fossil spectra all have generic 
features which can be understood from a physical stand- 
point. In Fig.O we show results both from volume weighting 
the individual spectra of SPH particles as well as the median 
CRe spectrum. The normalization of the two spectra differ 
significantly, where the volume weighted spectrum is dom- 
inated by about 1/5 of the particles which also make up a 
similar fraction of the total volume in the cluster outskirts, 
except for high momenta (p > 10^) where the particle frac- 
tion that dominates the volume weighted spectrum becomes 

^ For each simulated cluster, the CRe spectrum is followed on 500 
randomly selected SPH particles. Convergence studies show only 
marginal differences with increasing number of SPH particles. 
^ We chose a radial range of 1.3-1.5 iJ200 because it coincides 
with the average location of cluster accretion shocks, which are 
thought to light up reli cs (in fact, simulat ions show it can be even 
larger - e.g., Fig 10 of lVazza et"al] 1 I2OIII) .') 
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Figure 3. Representative electron momentum spectrum /cRc of 
a post-merging cluster in the virial region between (1.3—1.5) -R200. 
Top: we show the median injected spectrum (red dashed), the 
median cooled spectrum (dark blue solid) and its 68 percentiles 
(light blue dotted) . O verplotted are representative sample spectra 
of individual SPH particles (grey dash-dotted). The total volume- 
weighted spectrum (purple dash-dotted). The cosmic rays are 
compared to the thermal population described by a Maxwellian 
with a temperature of A;T2 = 3 keV. Bottom: spectral index of the 
median electron spectrum. 

progressively smaller. In this work, however, we take a con- 
servative approach, and focus on the median spectrum. Note 
that the cooled spectrum deviates significantly from the in- 
jected spectrum (shown with a red dashed line). The cooled 
spectrum can be characterized by four regimes: (1) the sub- 
relativistic Coulomb cooling regime (p < 1) where the CRs 
injected in the most recent shocks dominate the population; 
(2) the relativistic Coulomb cooling regime, which includes a 
peak at p ~ 1) due to the transition from non-relativistic to 
relativistic cooling; (3) the adiabatic regime (10 <p < 10^) 
where the cooling time is long (~ few Gyrs); (4) the radia- 
tive cooling regime (p > 10^) where the synchrotron and 
inverse Compton losses start to steepen the CRe spectrum. 

In the lower panel of Fig. |3] we show the spectral in- 
dex of the CRe spectrum, which is the same as that for 
the CRp. The injected spectra for CRp's has a concave 
shape with a spectral index ~ 2.8 at p ~ 1 which has fiat- 
tened to ~ 2.3 at p ~ 10^. This spectral shape is a conse- 
quence of the cosmological Mach n umber distribution tha t 
is mapped onto the CR spectrum jPfrommer et al.1l2006l ). 
It can be understood qualitatively as follows. The charac- 
teristic Mach number declines with time in a ACDM uni- 
verse, due to the slow-down in structure formatiorQ. Thus, 
shocks at early times have higher Mach numbers and harder 

^ By contrast, in an Einstei n-de Sitter Clm = 1 universe, struc- 
ture formation is self-similar llBertschi ngcr"l985). and the cosmic 
Mach number should not show any evolution, aside from non- 
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-1 



Ta = lO^K 
ing a = 2.5, a:inj,2 



spectra; late time shocks have softer spectra. However, early 
shocks also have lower normalizations, because they take 
place in a colder medium. Specifically, pinj oc T^''^— the 
power-law attaches to the thermal distribution at lower mo- 
menta in a colder medium. Thus, the normalization C oc 

C2/C1 « (10^)°-^5(3,.^.^2/2._^)i.5 ^ ^00, adopt- 
Xinj,!, and assuming a constant num- 
ber density of thermal electrons for this order-of-magnitude 
argument. We shall soon see ( tj3.2l and Fig.[S} that the com- 
peting effects of structure formation and the expansion of 
the universe conspire to keep gas density that ends up in 
the cluster outskirts relatively constant, so adiabatic evolu- 
tion does not significantly alter this conclusion. The upshot 
is that late-time weak shocks dominate at low p due to their 
higher normalization (hence, the spectral slope at low p is 
softer), while early-time strong shocks dominate at high p 
due to their harder spectra (hence, the slope here is harder). 
What is important to note is that the spectral index in the 
adiabatic regime p ~ 10 — 100 is relatively soft {a « 2.5), 
implying that primarily weak shocks M ~ 3 contribute. 

It is also instructive to consider the relative contribu- 
tion of different cosmological epochs to the spectra, shown 
in Fig. [4] There are two important effects: the decreasing 
normalization of spectra injected at early times (for the 
reasons discussed above), and the strong effects of cooling. 
Both conspire to decrease the importance of early time in- 
jection. Thus, the first effect means that most electrons in 
the peak and 'adiabatic' regime 1 p ^ 100 came from the 
last ~ 3 Gyr. In the regime p ^ 1, where sub-relativistic 
Coulomb cooling dominates, the contribution is even more 
recent-essentially within the last cooling time, ~ 0(10*) yrs. 
We now turn to understanding the shape of this spectrum 
in detail. 



3.2 Analytic Model for Fossil Electrons 

We develop an analytic model for the fossil distribution 
function /cro(p) in the presence of multiple cooling pro- 
cesses. There are two clear limiting cases: an "impulsive in- 
jection" scenario where an initial population of relativistic 
particles cools passively, and a "steady injection" scenario 
where a steady state balance between injection and cooling 
is achieved. The "impulsive" scenario is appropriate when 
the timescale between shocks is greater than the cooling 
time. The finite time resolution of our post-processed sim- 
ulations means that our simulation results are a linear su- 
perposition of such impulsive solutions. The "steady state" 
scenario is appropriate when the injection rate of relativistic 
particles is fairly constant, and tcooi(p) <C in (where in is 
the Ifubble time) so that the population equilibrates on a 
short timescale. 

Naively, since the CRs are injected at discrete shocks, 
one might expect the impulsive approximation to be most 
appropriate. In fact, the opposite is true. The timescale 
on which the shock injection process changes is the dy- 
namical time, tin] ~ tdyn ~ tu/y/A ~ IAA^q^^ Gyr, 
where tn ~ 10 Gyr, A50 = A/50 and A = p/p. On, the 




Figure 4. Diflerential spectral build up of the median electron 
distribution function, /cRe, as a function of cosmic time, for a 
post-merging cluster in the region between (1.3 — 1.5) R200- The 
black solid line shows the total /cRe- The thin grey lines show 
the individual differential contributions to /cRc from different 
times, where we highlight the contributions spaced by Gyr in- 
tervals with different colors. The spectral features are shaped by 
different cooling processes: non-relativistic Coulomb cooling for 
p < 1, a broad feature indicating the transition from relativistic 
Coulomb cooling to the adiabatic regime for 1 < p < lO'^, and 
inverse Compton cooling at higher momenta. 



other hand, as we shall soon see, the main regime of in- 
terest where cooling significantly modifies the spectrum is 
Coulomb cooling in the trans- and sub-relativistic regime, 
where icooi ~ lp^(no/2 x 10~^cm~^)~^ Gyr for p < 1. 
The fact that tcooi ~ iinj implies that the injection rate is 
roughly constant on the cooling timescales over which the 
population equilibrates. Moreover, the Coulomb cooling rate 
itself — which depends on the gas density — changes on the 
dynamical time tdyn. Figure [5] shows the Lagrangian evolu- 
tion of the physical gas density that ends up in the virial 
region of a cluster. The presence of clumping implies a non- 
Gaussian density distribution that biases the mean of the 
distribution upwards in comparison to the median, the lat- 
ter of which also shows a much smoother density evolution. 
In fact, the physical density does not even evolve signifi- 
cantly on timescales t 3> tdyijfl- Thus, on the short equilibra- 
tion timescale tcooi for the sub and trans- relativistic regimes, 
the injection and cooling rates are roughly constant, and a 
steady state solution is valid. However, before constructing 
these solutions, it is worthwhile to examine the injection and 
cooling processes in more detail. 

The right-hand panel of Fig. [5] shows the Lagrangian 
evolution of the median CRe injection rate of gas ending 
up in the outskirts, (1.3 — 1.5) R200, of cluster g72a, which 
experienced a merger 2 Gyrs ago. We can see that it is fairly 
steady over cosmological timescales. Similarly, we have also 
found that the injection rate of a cool core cluster that has 
not experienced a major merger in the last 7 Gyrs is either 
decreasing slowly or constant with time. Note that this is 



gravitational events such as reionizati on and effec ts from finite 
mass resolution in the simulations (iPfrommer et al.ii2006,) . 



For want of a better explanation, we view this as coincidental 
cancellation between the competing effects of Hubble expansion 
and structure formation. 
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Figure 5. Lagrangian evolution of the thermal electron number density, rtc (left), and the injection rate density of CR electrons, Qinj^ad 
(right). We consider a representative sample of SPH particles that end up at the virial region of a massive, post-merging cluster (g72a) 
between (1.3 — 1.5) -R200 at ^ = and follow their evolution as a function of time since the Big Bang. Left: the mean (blue dotted) 
and median (red solid) ric. Their difference is due to its non-Gaussian distribution caused by clumping, which biases the mean upwards. 
Right: we show Qinj,ad fof three different momenta: p = 0.1 (green dotted), p = 1 (red solid, weighted by factor 10^), and p = 100 (blue 
dashed, weighted by factor 10^). The three small insets in each panel show CR electrons with a momentum p > 1 at different times; 
the density distribution of the fossil CR electrons (left), and the injection rate density (right). The variance is caused by different shock 
strengths, injection efficiencies, thermal densities and decreases with time. 



only true for the median spectrum (or any other spectrum 
which is averaged over a large number of fluid elements). An 
individual fluid element can have a more stochastic injection 
history, and the fossil electron population in the outskirts of 
the cluster can vary spatially. This implies that when a weak 
shock propagates across a cluster, not all regions will hght 
up with equal intensity. This spread in the CRe population 
is shown in the insets on the left panel of Fig. [5] while the 
spread in the injection rate is shown in the right panel. These 
distributions narrow with time, but can still span 1-2 dex at 
z = 0. 

For now, we focus on understanding the median spec- 
trum. Figure [S] shows that the median density over the last 
~ Gyr is rio ~ 1 X 10~^ cm~^ (corresponding to A ~ 100), 
while the median CRe injection rate density is ncRe(p > 

1) = /rQinj,ad(p)dp ~ Qinj,ad(p = 1) ~ IQ-'^Cm-^S-V 

We can understand this from simple order of magnitude 
arguments. From Fig. [3l the spectral index for the adi- 
abatic spectrum is a ~ 2.5, corres ponding to .M ~ 3 
shocks. This is consistent with Fig. 1 of lKang fc Rvul |20l3), 
where they find that the kinetic energy flux of shocks has 
a sharp drop-off for shocks with A4 Si 3. A M 3 shock 
injects CRs at the cluster periphery with number den- 
sity ncRe(p > 1) ~ ne??CR,lm(p*/Pinj)~" + ^ ~ 10"" Cm"^, 

where 7?cR,iin = ncn/rie ~ 10~^ is the injection efficiency 
for = 3 shock (where pinj ~ 4ptherm), and the factor 
of (p*/pinj)~"^^ converts n(p > pinj) to n(p > p*), and 
p, ~ 1. Since the shocks operate on a dynamical timescale 
tdyn ~ fn/x/A ~ 1 Gyr (A/100)~^^^, the estimated injection 
rate is ncR(p > 1) ~ ncR,o(p > l)/^dyn ~ 10"^^ cm~^ s~^, 
which agrees with the simulation results quite nicely. 

As for cooling. Fig. [1] illustrates the cooling time as a 
function of energy, for ric ~ 10~^cm~"^ at 2 = 0. There is a 
quasi-adiabatic regime pi < p < P2 where pi ~ 10, P2 ~ 10^ 
where tcooi > in; otherwise, for p < pi, Coulomb cooling 



dominates, while for p > p2, inverse Compton cooling domi- 
nates. For our purposes, Coulomb cooling is the most im- 
portant cooling process. Inverse Compton cooling merely 
shifts particles from the high energy tail to the adiabatic 
regime (which acts as a 'road-block'), where they are still 
available for reacceleration. Moreover, the relative number 
of affected particles in the power-law tail is relatively small. 
By contrast. Coulomb cooling affects the low-energy regime 
where most of the particles are by number, and can cause 
them to migrate to low thermal momenta, where they are no 
longer available for reacceleration. Note that the Coulomb 
cooling function (equation I24|) undergoes a rapid transition 
from relativistic energies (6c ~ const) to sub-relativistic en- 
ergies (6c oc p~^), which implies cooling times tcooi oc p (for 
1 < p < 30, where the upper limit depends on density and 
redshift) and tcooi oc p^ (for p < 1). This rapid thermaliza- 
tion of sub-relativistic particles creates a peak in the steady 
state spectrum at p ~ 1. 

Let us now consider a toy model where the injection 
and Coulomb cooling rates are constant. As previously men- 
tioned, this is a reasonable approximation, as they both 
change on a dynamical time tdyn which is longer than the 
equilibration time tcooi. For the same reason, we can ig- 
nore adiabatic changes to the electron spectrum. The mag- 
netic fields which govern synchrotron cooling also change 
on a dynamical time (for adiabatic changes, or a turbu- 
lent dynamo), while inverse Compton cooling changes on 
cosmological timescales. It is straightforward to solve for 
the steady state solution, using the Vlasov equation. How- 
ever, this is only appropriate in the low and high energy 
regimes, where cooling (and hence equilibration) timescales 
are short. In the adiabatic regime (pi < p < P2), the elec- 
tron population is time-dependent — it simply grows with 
time. Fortunately, when the injection and cooling function 
are time-independent, one can construct self-similar solu- 
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Figure 6. Comparison of the self-similar model (red dotted) of 
equation II28I I to the median CRe spectra, /cRoi oi our simulated 
post-merging cluster g72a. We show the median /cRo for two 
different time resolutions, 100 Myr (dark blue solid) and 1 Gyr 
(light blue dashed), indicating impressive convergence at p > 100, 
which is the relevant energy regime for radio relic emission. 



tions which co nnect the steady and non-steady populations 
llSarazinlligggl ). In these solutions, the overall normalization 
varies with time, but the shape of the energy spectrum re- 
mains the same if p is scaled by some characteristic value. 
It is: 



/sclf-sim(p, t) 
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(28) 



The approximation in equation H28p is valid as long as pic 3> 
PCoul, where pic = l/(6ic,o At), and pcoui = 2bc(p,t) At|f|. 
This expression explicitly makes the approximation that the 
injection rate Qinj.ad and cooling rates be, he are indepen- 
dent of time. It thus only requires two time-independent 
input parameters, Qinj,ad,nc. In Appendix [Cl we provide a 
power- law model for Qinj,ad with the spectral index ainj. In 
Fig. |6] we compare the results of our simulations of g72a 
(which have 100 Myr time resolution) to the self-similar 
model in equation (|28p . The agreement is remarkable, given 
that injection and cooling rates are not exactly constant. We 
also show a calculation with lower time-resolution At = 1 
Gyr, where the agreement is worse, especially in the mo- 
mentum regime where IC cooling is substantial. The rea- 
son for this is that there are less SPH particles, which got 
recently injected. Hence IC cooling modified the spectrum 
more severely and removed all CRes with a momentum 



^ We adopt a characteristic timescale for the cooling processes 
of At ~ 9 Gyrs (i.e. pic ^ 30 and pcoul ~ 4.5), which roughly 
corresponds to the formation time of a galaxy cluster in our sim- 
ulations. For timescales greater than a few Gyrs, we introduce a 
factor 2 in pcoul to effectively account for the order unity decrease 
in density in the cluster outskirts. 



Figure 7. We show the contribution per logarithmic interval in 
momentum to the integral of the downstream distribution func- 
tion of reaccelerated CR electrons in equation II30I I. For a strong 
shock {A4 > 4), there is a broad peak of fossil electrons with 
1 < p < 100 dominate the contribution to reaccelerated CR elec- 
tron distribution, while for a weak shock, higher momenta with 
p > 10 contribute equally. 



above pic- As a result, median spectrum is biased low at 
high energies. We shall return to this issue of finite time 
resolution in M.2\ 



4 REACCELERATION 

4.1 Which Momentum Regime Matters? 

Downstream of a shock, t he CR distrib ution function can 
be written as (Bell 1978bl: lDrur"^ll983bh : 

rp 

/2(p) = (a + 2)p'°' / p'""Vo(p')dp' 



+ — ) H{p-pini), 



(29) 



where a is the test particle power-law slope and H{p — pinj) 
is the Heaviside step function. The first and second terms 
on the right hand side refer to the contribution of relic and 
freshly injected CRs respectively; we focus on understanding 
the former in this section. This first t erm can be w ritten in 
a more physically instructive way as (|Drurvlll983bl ): 

/r^^"'(p) = + 2) (^) " .Mp')H{p - p')d(logp') . 

(30) 

Thus, the downstream spectrum is a convolution of the up- 
stream spectrum with a truncated power-law. Note that 
since particles are conserved and do not lose energy dur- 
ing the acceleration process, we expect 712'^'"^'^° (p > pinj) — 
rno'^^^^ip > Pinj)i where r = P2/P0 is the shock compres- 
sion ratio. Indeed one can derive equation (|30p from this 
assumption. These statements are reminiscent of the man- 
ner in which a thermal Maxwellian tail with p > pinj is 
converted into a power-law — nothing about the acceleration 
process relies upon a Maxwellian distribution, and the same 
conclusions hold for a non-Maxwellian tail. 
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For concreteness, an illustrative example can be use- 
ful. Consider the case when the relic s pectrum is a power- 
law fo(p) = h. ,{ph,)-'^H{p - p,) l|Kang fc Rvufl201ll : 
iKang et ahl I2OI2I'). Substituting fo jp) into equation ((30)) 



yields for a 7^ /3 iKang fc Rvull201ll ): 



/2(P) 



Q + 2 



/o(p) 



(31) 



(32) 



where 7 = min(a,/3). Thus, for p ^ p* , the reaccelerated 
distribution function asymptotically becomes a power-law, 
with a spectral slope corresponding to the shallower (i.e., 
harder) of the initial spectrum and the reaccelerating shock. 
Of course, our fossil electron spectrum is very different from 
a power-law, due the effects of cooling: it is a sharply peaked 
function, with a continuously varying slope. Given that our 
simulations are only accurate in a limited momentum range 
p', we can ask: what range of p' must crucially be resolved? 
Naively, one might assume that that since reacceleration 
conserves number density, it suffices to resolve the peak of 
p'/(p'), where most particles reside. However, the additional 
weighting by a power-law in equation pO|> means that this 
condition is insufficient: although p'/(p') robustly peaks at 
p' ~ 1, electrons here may be 'too far from the action', once 
multiplied by the lever arm (pcmit/p')~"- 

Let us consider which initial momentum regime con- 
tributes most to observed synchrotron emission of reaccel- 
erated fossils. The characteristic synchrotron freauenc-J"1 is 
~ 37^ fc, where fc is the non-relativistic cyclotron frequency; 
this implies that for a given observation frequency v^, the 
greatest contribution comes from electrons with 



Pcmit ^ O'cmit 



5 X 10'' 



IGHz 



1/2 



B 
5^G 



-1/2 



(1 + 2 



,1/2 



(33) 

where for instance B « 5/iG has been inferred from high 
resol ution measurements of spectral aging in the sausage 
relic (Ivan Weeren et al]|2010l ). Figure [7] illustrates the con- 
tribution per logarithmic interval to the integral in equa- 
tion pO|) . for Pcmit = 10*. For a strong shock, we see that 
the integrand has a broad peak ranging from 1 ;S p' 100, 
while for a weak shock it increases monotonically with p', 
receiving its dominant contribution from 10 



2 < „' < 



10* 



We can understand this as follows. The integrand peaks 
when a — /3(p') = 0, where /?(p') = — dlog /o/dlog p', i.e. 
when the power-law slope of reaccelerated particles and the 
fossil distribution function have the same slope. It picks up 
most of its contributions from a neighborhood of this re- 
gion. From Fig. |3l we see that /?(p') ranges from zero (due 
to NR Coulomb cooling) to extremely steep positive slopes 
/3(p') > 3 (due to inverse Compton and synchrotron cool- 



1'^ Here we use the mon ochromatic approximatio n of synchrotron 
emission (see App. B of lEnfilin &: Sunvaevll2002l ) where the syn- 
chrotron kernel is replaced by a delta distribution, S(v — Vs) with 
Va = 3e_B sin S 7^ / (27rme c) , where 9 is the CRes' pitch angle and 
which gives the exact synchrotron formula for a power-law elec- 
tron population with spectral index a = 3 and attains only or- 
der unity corrections (< 20 per cent) for small spectral changes 
Aa < 0.5. 



ing), this tangent is guaranteed to exist. In practice, the 
relevant regime where 2 < /3(p') < 3 lies in 1 ^ p' ^ 10'', 
the 'adiabatic' regime which only only mildly affected by 
relativistic Coulomb cooling or inverse Compton cooling. 
Regions which have much shorter cooling times and con- 
sequently had their injected slopes strongly modified do not 
contribute significantly. This suggests the fortunate con- 
clusion that we do not need to accurately resolve regions 
strongly affected by cooling to accurately predict the reaccel- 
erated spectrum. An important point to note is that in this 
crucial adiabatic momentum regime 1 ^ p ^ 100, the spec- 
tral index of the injected spectrum a « 2.5 (e.g., see bottom 
panel of Fig. [3]), indicating that weak shocks {A4 ~ 3) dom- 
inate the assembly of the fossil CRe spectrum here. We will 
return to this point in H5.ll 



4.2 Do Our Simulations Have Sufficient Time 
Resolution? 

We now turn directly to the question of required time res- 
olution. Our finding in H3.2I that injection and cooling is 
roughly time steady appears to contradict the approach we 
have taken in our simulations, where we inject CRes in dis- 
crete bursts (due to finite time resolution). In Fig. |S] we 
show the distribution function obtained by mimicking the 
same procedure we followed in our simulations: injecting 
the electrons in discrete bursts, and then passively cool- 
ing the injected population from each time step, via equa- 
tions (IA13|) - (IA15|) . for different values of At (in our simula- 
tions. At = (lOOMyrs - 1.5Gyrs), see Table[T]for details). 
For a fair comparison, we have assumed that cooling and 
injection are exactly constant. By comparison with the left 
panel in Fig. [Sj we see that in regions where At < tcooi, our 
computational procedure correctly approaches the analytic 
solution. However, when At > tcooi, the discrete grid over- 
estimates the importance of cooling, causing the numerical 
solution to fall below the analytic one. 

However, our results from the previous section suggest 
that the regions where cooling is important do not con- 
tribute significantly to /(pcmit). In particular, the main im- 
portant mechanism is non-relativistic Coulomb cooling at 
p ^ 1. This momentum regime of the fossil spectrum does 
not contribute significantly to /(pcmit) (as in Fig.[7| because 
it is 'too far from the action' at Pcmit, and also has been 
significantly depopulated. Thus, our finite time resolution 
does not significantly affect results. We show this explicitly 
in the right panel in Fig. (8] which shows pcmit ./'rcacc(pcmit) 
from the reaccelerated population for the various curves in 
the left panel of Fig.O as a function of shock Mach number. 
The differences are small, given other uncertainties in the 
problem. This shows that the contribution from fossil reac- 
celeration can indeed be estimated accurately analytically. 
Simulations are nonetheless invaluable for their ability to 
shed light on spatial (see the insets of Fig. [5} and tempo- 
ral variations (see ^5.'2\i in injection and cooling rates, all 
of which can lead to significant scatter in the CRe fossil 
population. 



© 2011 RAS, MNRAS 000.[Tlf23l 



12 A. Pinzke, S.P. Oh and C. Pfrommer 




W- 10"' 10° 10' 10- 10- 10* 2 3 4 5 6 7 8 

p = PJ mj: M 



Figure 8. Continuous versus discrete injection. We sliow the cooled fossil CR electron distribution function as a function of 
normalized momentum. The self similar distribution (red dotted) assumes constant continuous injection and constant density and 
is compared to discrete injection scenarios of CRes with different duty cycles: 1 Myr (dark orange), 10 Myr (light orange), 100 Myr 
(light green), 1 Gyr (green), 5 Gyr (dark green). We assume an injected CR electron population given by equation IjClj l and an 
electron number density, rie = 10~^cm~''. The distribution function from discrete injection approaches the self similar solution 
when the injection timescale becomes smaller than the cooling time. Note that for this figure we adopt a characteristic assembly 
time of the CRes of 10 Gyrs, hence the 5 Gyrs curve (at p = 10) and the 1 Gyr curve (at p = 10'') are superpositions of two and 
ten curves, respectively, which explains the visible shoulders for these curves. 
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Figure 9. Median p/cRc for direct injection (dashed line) and 
reacceleration (solid line), for different Mach numbers (color 
coded); also shown is the fossil distribution function (dot- 
dashed). /cRc for direct injection is a power-law attached to 
the Maxweilian, while /cRc- The distribution function falls much 
more drastically at low Mach number for direct injection (dashed) 
than for reacceleration (solid). 



4.3 Principal Result: Reacceleration Dominates 
over Direct Injection 

Thus far, we have focused on accurate calculations of the 
fossil electron distribution function. We now explicitly illus- 



trate, using equation (|3Up . how this population of electrons 
is transformed by a shock. Figure [9] shows the initial fossil 
electron distribution function and the transformed distribu- 
tion function after Mach number A4 = 2, 3, 8 shocks; the 
population due to direct injection is also shown. The distri- 
bution function for direct injection is a power-law attached 
to the Maxweilian, while the reaccelerated distribution func- 
tion is a boosted version of the fossil distribution function, 
with an asymptotic power-law tail. This tail has the same 
power-law slope a as for direct injection. Note that while 
the total number of fossil electrons is conserved, the number 
density of reaccelerated electrons increases over that of the 
fossil population by the shock compression factor r, which 
explains the normalization boost upon reacceleration. It is 
clear that while the electrons from direct injection have a 
strong Mach number dependence, this trend is much weaker 
for the reaccelerated population. This suggests that reaccel- 
eration could dominate at low Mach numbers. 

Figure [TO] shows the main result of this paper, the ra- 
tio /rcacc(Pcmit)/,finj(Pemit), Or the ratio of the distribution 
function at the primary emitting frequency for the reac- 
celerated population to the same for the freshly injected 
population (dark blue curve). We see that at high Mach 
numbers {A4 > 4), relic reacceleration and fresh injection 
are comparable, but for low Mach numbers {M < 4), fos- 
sil reacceleration vastly dominates over fresh injection. We 
also compare the reacceleration and direct injection models 
to a straw man model ('M-const') where a constant frac- 
tion of the shock energy Cinj.c ~ 0.01 is converted to rela- 
tivistic electrons, independent of Mach number. As previ- 
ously mentioned, while we do not believe this assumption 
to be physically realistic, it has been widely used in the lit- 
erature, and gives us a baseline for comparison. While the 
'M-const' model gives comparable results to the other two 
models for > 4, it differs sharply from the direct injec- 
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tion model for A4 < 4, over-predicting the CRe population. 
This makes sense, since physically we expect the injection 
efficiency to plummet at low Mach numbers. By contrast, 
the 'M-const' model is similar to the reacceleration model 
down to yVf ~ 3, but for M < 3 it under-predicts the CRe 
population. This has two important consequences. Firstly, 
it means that 'M-const' models have been getting the right 
answer for the wrong reasons: while they can match observa- 
tions of ^ 3 relics, this is not necessarily because acceler- 
ation efficiency is independent of Mach number. Instead, our 
model provides a physical basis for the observed brightness of 
low Mach number relics, through the existence of fossil CRe 
(see also Fig. llSp . Secondly, at extremely low Mach numbers 
< 3, reacceleration predicts many more CRe than even 
this model. Thus, we predict many more relics with steep 
spectra which are potentially observable with LOFAR (see 

How can we understand the Mach number dependence 
of yrcacc (Pemit )//inj(Pcmit)? It is casy to understand why it 
rises steeply at low Mach numbers. The Mach number de- 
pendence of /rcacc(Pcmit) is solcly duc to the dependence of 
the power-law slope a on Mach number; the number of par- 
ticles available for acceleration is fixed. As a steepens at 
low Mach number, /reacc(Pcmit) falls in a power-law fashion. 
This mostly happens for A4 ^ 4, when a begins to evolve 
significantly. On the other hand, fresh injection is affected 
both by the change in a and more importantly, the change in 
the number of particles available for acceleration. As can be 
BEisily derived from the jump conditions, low Mach number 
shocks have lower thermalization efficiencies: the post-shock 
gas has a lower temperature and higher bulk velocity. This 
means that Xinj = pinj /pthcrm has to be higher for a particle 
to be able to overcome the bulk fluid motion to cross the 
shock. This increase in Xinj on the Maxwellian tail exponen- 
tially reduces the number of particles which are accelerated 
at low Mach number shocks, which is why /inj(Pcmit) is ex- 
ponentially suppressed. 

On the other hand, the fact that /rcacc(Pcmit) and 
/inj(Pcmit) are roughly comparable at high Mach numbers 
may appear somewhat surprising. Since the build-up of the 
CRe fossil population is due to the interaction of continuous 
injection via multiple shocks with cooling processes, a priori 
we might expect that it could potentially be orders of mag- 
nitude smaller or larger than that due to direct injection at 
a single shock. In particular, as we have seen in Fig. [T] since 
much of the contribution to /reacc(Pemit) comes from an adi- 
abatic regime which grows monotonically on cosmological 
timescales, we might expect the reservoir of fossil CRe built 
up to outweigh the contribution from a single shock. In fact, 
as seen in Fig. |3] and Fig. [S] the injection rate is fairly time- 
steady, with late times having a somewhat larger contribu- 
tion since the gas is getting hotter, increasing pthcrm. The 
characteristic injection timescale is t^yn ~ tn/y/A ~ 1 Gyr, 
implying that the fossil CR electron population should be 
tH/tdyn ~ 10 times larger. The fact that cooling operates, 
and the reduced contribution from early times makes the 
value of /reacc(Pemit)//inj(Pemit) ~ 3 showu lu Fig. [lO] reason- 
able. We also note that the fossil CRe population can show 
significant scatter from cluster to cluster, up to an order of 
magnitude (which we shall discuss in !j5]2} , so one should not 
over-interpret the exact value of /rcacc(Pcmit)//inj(Pemit). 
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M 

Figure 10. The ratio of the distribution function evaluated 
at p = Po/moC = 10^ (dotted) and 10'' (solid), for differ- 
ent acceleration models. The fiducial dark blue curves show 
/cR(reacc.)//cR(inj). Reacceleration is roughly comparable to 
direct injection at high Mach numbers and vastly exceeds direct 
injection at low Mach numbers. The green (red) curves com- 
pare reacceleration (direct injection) with a strawman model- 
commonly used in the literature-where a constant fraction 
Cinj.e = 0.01 of the shock energy goes into electrons, indepen- 
dent of Mach number. The straw man model has roughly similar 
behavior to the other models at high Mach number, but vastly 
divergent behavior at low Mach number. Note that the reaccel- 
erated CRe fossil population, /cR(reacc.), is derived from the 
self-similar model in equation II28I I using a characteristic cooling 
time At ~ 9 Gyrs. 



5 VARIATIONS IN THE FOSSIL ELECTRON 
SPECTRUM 

5.1 Dependence on Shock Acceleration Model 

In ij2.2.3l we discussed 3 different models for particle ac- 
celeration, but thus far largely used one model ('M-aR') to 
compute the fossil electron spectrum. Fig. [10] has already 
highlighted differences between the 'M-aR' model and the 
'M-const' model for a single shock. We now propagate these 
differences between the three models when the entire assem- 
bly history of the fossil population is taken into account. 

The resulting fossil electron population for the test par- 
ticle model ('M-testp'), the CRe spectral index renormaliza- 
tion model ('M-aR'), and the model with constant acceler- 
ation efficiency Cinj,c = 1 per cent ('M-const') is shown in 
Fig. [Til We see that while the 'M-aR' and 'M-testp' models 
give very similar results, the 'M-const' model yields a fossil 
distribution function with a significantly higher normaliza- 
tion. 

These results are easy to understand. From the bottom 
panel of Fig. [3] (and discussed in M.ip . we see that the spec- 
tral index of the median injected population before cooling is 
Oe ~ 2.5, i.e., weak shocks with ~ 3 dominate the assem- 
bly of the CRe fossil population. The 'M-aR and 'M-testp' 
models are identical at low Mach numbers; they only differ in 
how energy conservation is enforced at high Mach numbers 
(even so, differences are slight). By contrast, the 'M-const' 
model is broadly similar to the other models at high Mach 
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Figure 11. Comparison of median CRe spectra in the cluster 
outskirts of g72a at z = in different acceleration models. Shown 
are the test particle model ('M-testp'), the CRe spectral index 
renormalization model ('M-aR'), and the model with constant 
acceleration efficiency Cinj,o = O.Oi ('M-const'). 

numbers, but diverges sharply at low Mach numbers, with 
a significantly higher normalization (c.f. the red curve in 
Fig, nop , since by construction Cc.inj is constant, whereas it 
plummets drastically at low A4 for the other models. Since 
fossil CRes are primarily put in place by weak shocks, it is 
not surprising that they are much more abundant in this 
model. This amounts then to an increase in normalization 
of the 'M-const' model by roughly an order of magnitude in 
comparison to the other models (as can be inferred for the 
acceleration efficiency at = 3 shocks from Figs. 1101 and [2l 
considering that the relevant shocks are preferentially those 
at late time in a hot medium with particle energies of a few 
keV). 

These results indicate our main conclusion that fossil 
reacceleration must dominate at low Mach number is inde- 
pendent of the shock acceleration model. If one attempts to 
explain the observed brightness of weak shocks by modifying 
shock physics and increasing the acceleration efficiencies at 
low M, the same would hold for previous weak shocks and 
the normalization of the fossil CRe population would be cor- 
respondingly larger. Since the shock acceleration efficiency 
largely cancels out, the relative importance of reacceleration 
vs. direct injection depends on accretion/merger history and 
cooling, for which (unlike shock acceleration) there is rela- 
tively little uncertainty in the physics. 

5.2 Cluster-to-Cluster variations 

7s the build-up of a CRe spectrum m cluster outskirts a uni- 
versal process? To answer this question, we compare the CRe 
distribution for three typical clusters at different times in 
Fig. 1121 a small cool core cluster, a large cool core cluster, 
and a large post-merger cluster. It is remarkable how similar 
the spectral shape is between the different clusters and at 
different times. This generic shape can be understood from 
the analytic model presented in H3.2I The peak is at p ~ 1 
and the normalization for each cluster changes by a factor 
less than 10. The recent merger cluster g72a shows the high- 



est CRe number density, while the CRes in the two cool core 
clusters have a lower abundance. 

What is the origin of this dispersion? Does it originate 
from differences in CRe injection or cooling, between clus- 
ters? In Fig. 1131 we show the CRe spectra for all 14 clusters, 
both with (blue) and without (green) cooling. The adiabatic 
spectra (which were only transported adiabatically with- 
out accounting for cooling) exhibit remarkably little scat- 
ter, while the cooled spectra show considerably more scat- 
ter. This suggests that the variance between clusters arises 
from the cooling rather than the injection process (for more 
on scatter in the adiabatic spectrum, see Appendix [Cl and 
Fig. [en . 

Two other features in Fig. [13] are worth noting. Firstly, 
cooling introduces increased scatter even in the ostensibly 
quasi-adiabatic regime 10 ^ p ^ 100. Note that this regime 
has fcooi > tn only for Ue ^ 10~* cm~^ and at 2 = 0. For 
higher gas densities and redshifts, some cooling is possible; 
thus, variations in gas clumping (which affects the amount 
of Coulomb cooling) and shock history (which affects the to- 
tal amount of cooling; cf. the difference between "impulsive" 
and "continuous injection" solutions) . The fact that there is 
substantial scatter at p ~ 100 (where only inverse Comp- 
ton cooling is important) points toward the latter. Namely, 
clusters which receive their dose of CRes earlier than others 
undergo more inverse Compton cooling. The different build- 
up and normalization of the CRe population in the cool core 
and recent merger clusters in Fig. [T^] also supports this. Sec- 
ondly, the scatter at low p is considerably higher than at high 
p. Part of this is because the cooling time is shorter at low 
p, and hence differences in shock history-which are averaged 
over shorter timescales-are amplified. However, part of the 
difference is artificial, and due to our finite time resolution, 
which overestimates the effects of cooling ( H4.2l and Fig. [S}. 
For instance, the two clusters with the higher normalizations 
in the Coulomb cooling regime are g72a and g72b, which 
have a time resolution ten times better than the other clus- 
ters. Fortunately, as shown in ii4.1l this momentum regime 
is unimportant in its contribution to the reaccelerated spec- 
trum at Pemit. By contrast, the cooling time is generally 
longer than our time resolution in the p ^ 10 regime. 

Finally, by integrating over momentum f equation I30|) . 
we show the scatter in the observationally more relevant 
quantity pcmit/(pomit), where pemit ~ lO", in Fig.[l4l plotted 
against cluster mass M200. The scatter is about an order 
of magnitude, and somewhat below the apparent scatter in 
Fig. 1131 since it is weighted towards the lower scatter, high 
p regime. There may be a trend for the CRe distribution 
function to increase with M200 for merging clusters, though 
our sample size is too small to make definitive statements. 



6 COMPARISON WITH RADIO RELIC 
OBSERVATIONS 

In this section we derive the radio synchrotron emission from 
our fiducial 'M-aR' CRe model and compare it to observa- 
tions of radio relics. These comparisons assume the relics 
lie between (1.3-1.5)J?2oo and explore a limited range of pa- 
rameter space; they are only meant to be illustrative. Our 
main task is to explore distinctive observational signatures 
of fossil electrons. As expected, low Mach number shocks are 
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Figure 12. Lagrangian evolution of CRe spectra of three different clusters. We show the median CRe distribution function of a 
representative sample of SPH particles that all end up at 2 = in the virial regions within (1.3 — 1.5) iJ200 (different time intervals 
since Big Bang are color coded). Shown are a large cool core cluster (gla, left), a cluster with a recent merger (g72a, middle), and 
a small cool core cluster (g676, right). There is strong evolution of /cRe in the post-merger cluster, which is noticeably reduced 
in the cool core clusters. While the normalization is increasing in our post-merger cluster that experienced a more violent recent 
history, the opposite trend is visible in the cool core clusters where CRe cooling processes dominate. However, the overall shape 
of the distribution function is very similar between those clusters that vary widely in dynamical stage and mass. 



Table 2. Relic sample. For simplicity we assume a downstream magnetic field of 5 fiG and a downstream temperature of kT2 = 
5 keV. 



relic name 


redshift 






thickness'"^) 


(4) 


P(1.4 GHz) 


reference 








[Mpc2] 


[kpc] 


[el cm""^] 


[10^2 erg s-l Hz-i] 




A 2256 


0.0594 


2.6 


0.6 


70 


6 X 10""' 


0.4 


Clarke et al. (2011) 


A 3667 


0.055 


4.7 


2.0 


62 


10-* 


4.1 


Rottgering et al. (1997) 


Sausage 


0.1921 


4.5 


1.5 


48 


3 X 10-5 


1.4 


van Weeren et al. (20101 


Toothbrush 


0.225 


4.6 


3.5 


45 


10-* 


6.0 


van Weeren et al. (2012) 


A 2744 


0.3080 


2.4 


2.6 


44 


4 X 10-* 


0.5 


Orrii et al. (2007) 



Notes: 

(1) The Mach number M is derived using the observed spectral index of the radio emission closest to the shock front. (2) Shock 
area, estimated from the largest linear size of relic squared. Note that there is a large uncertainty of the relic in the direction of 
the line of sight. (3) Relic thickness as estimated from the cooling length, equation I I35II . (4) Assumed electron density adjusted to 
match observed radio luminosities. 



much brighter if fossil electrons abound; we therefore pre- 
dict many more steep spectrum sources to be detectable at 
low flux limits. On the other hand, the relic luminosity func- 
tion is not a robust discriminant of models with and without 
fossil electrons. Thus, spectral information is needed to test 
our model. 

Given the many uncertainties, we adopt a simple model 
to estimate radio luminosities, which still takes the effects 
of cooling into account. The radio syn chrotron emissivity of 
a pow er-law distribution of CRes is l|Rvbicki fc LightmanI 
[l97i): 



12 



12 J 



where Jo 



32 Uc 
c (a + 1) 



(34) 



Here the radio spectral index — {a — l)/2, a is the spec- 
tral index of the CRe population, F is the gamma func- 
tion, and Vc ~ eB / {2nmec) is the cyclotron frequency. The 
~ 10 GeV electrons which emit at ~GHz frequencies (equa- 
tion [33| cool via IC and synchrotron emission over a post- 



shock distance: 

^cool ~ W2^cool 



200 



kpc 
B2 



U2 



5/iG 



103 kms-i 
1/2 / 02 

' i>cft,2 



IGHz 

-1 



-1/2 



(1 + ^) 



-1/2 



(35) 

where 2 = B| + Bomb, Bomb = 3.24(1 + zf y,G, and 
typical downstream velocities, tt2, in a relic 10^-10^ kms-\ 
The specific luminosity is given by: 

P. ~ (/cooi^s) X J{v) cc v"-!^ S^'^JI (36) 

^ + ^CMB 

where ^Ps is the shock area, and we assume that all post- 
shock variables are roughly constant over a distance /cool- 
For a given electron population, this simple estimate gives 
similar scalings to more ca reful calculations which inte- 
grate over the cooling layer ijHoeft fc Briiggenll2007l ). from 
the freshly accelerated to oldest electrons. In particular, 
we obtain the same result that cooling steepens the spec- 
tral index from (a — l)/2 to q/2. For weak B-fields 
<^ BcMB ~ 3.24(1 -I- z)^ fiG, the luminosity increases 
with B-field, P„ oc B^+"^^. However, this increase starts to 
saturate when B ^ Bcmb; for B^ > Bcmb, Pi' oc _B"/^-^ 
which is very weak (for instance, P^ oc B^'''* for a = 2.5). 
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Figure 13. CRe spectra in tlie virial regions of clusters at 
2 = 0. Shown are tlie median CRe spectra of all simulated clus- 
ters normalized to the same thermal electron number density 
'^c,norm = 10~^ cm"''; with cooling (light blue) and the median of 
those (dark blue), injected spectra (light green) and the median of 
those (dark green). The injected spectrum appears quite univer- 
sal with little scatter. By contrast, although the cooled spectra all 
have similar shapes, their normalizations can differ substantially, 
by up to a factor of ~ 50. 
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Figure 14. The momentum integrated quantity Pcmit/CPomit) 
(where Pemit ~ 10*) plotted against cluster mass. For each cluster 
the CRes are normalized to the same thermal electron number 
density no, norm = 10~^ cm~'^. The red crosses and blue X show 
the simulated cool core (CC) and merging clusters, respectively. 



This makes physical sense: when synchrotron emission dom- 
inates cooling, then all of the energy in relativistic electrons 
is emitted at radio wavelengths, independent of the B-field. 

Is the brightness of observed relics consistent with our 
calculated fossil electron population? This question is most 
accurately answered with the handful of observations with 
high spatial resolution or favorable geometry where the ef- 
fects of spectral aging can be resolved or otherwise mini- 
mized. We list these in Table [2l The spectral index closest 
to the shock front, = (a — l)/2, before cooling steep- 
ens the spectrum, is an accurate measure of shock Mach 



number. The shock area is estimated as ^'s ~ , where 
L is largest observed linear size of relic squared. For sim- 
plicity, we assume B2 ~ 5fiG (as for instance estimated 
for the Sausage relic using spatially reso lved observations of 
spectral ageing Ivan Weeren et aLll2010l ). and fcTa ~ 5keV, 
and only allow variations in Ue- Assuming the median fossil 
electron spectrum from our simulations, we can match the 
observed relic luminosity by reasonable variation^^ in the 
electron gas density Ue ~ 10~*cm~''. Thus, reaccelerated 
fossil electrons can clearly produce radio relics of the right 
luminosities. 

How do the different models compare in their predic- 
tions for radio luminosity as a function of Mach number? In 
Fig. 1151 we contrast the reacceleration model with the direct 
injection model ('M-aR'), where we adopt fiducial param- 
eters for shock area, temperature, density and B-field, but 
also take into account their possible spread (indicated by 
the shaded regions). We also show observations where the 
intrinsic spectral index can be accurately determined, as in 
Table [2] (filled diamonds; labelled; here ai, = {a — l)/2), 
and those where we correct for the effects of cooling (open 
diamonds; upper limits; here a^, = a/2, d e rived from the 
compilation of radio relics in iFeretti et all (|2012l )). While 
discrepancies with the fiducial curve in the reacceleration 
model (up to a factor of ~ 10 — 30) lie within the uncer- 
tainties, the much larger 100 — 1000^ discrepancies with 
the direct injection model means that reconciliation is wholly 
untenable. The radio flux from reaccelerated fossil electrons 
declines smoothly as a function of Mach number (factor 30 
difference between M ~ 2.5 and M ~ 4.5) compared to 
the direct injection scenario, where the flux is exponentially 
suppressed as the shock become weaker than ^ 4. This 
means that as more sensitive radio experiments are built, the 
reacceleration model predicts that there will be a substantial 
increase in the number of steep spectrum relics, especially 
since low M shocks are believed to be more abundant that 
the high M shocks. These results for the overwhelming dom- 
inance of reaccelerated fossils at low Mach number mirrors 
that in Fig. IIUI as it should. 

We also show on the same plot results for the constant 
acceleration efficiency = 0.01 ('M-const') model. Note 
how it does not suffer from exponential suppression at low 
Mach numbers by fiat, and to some extent mimics the be- 
havior of the reaccelerated fossil model. This ability to more 
or less fit the observations is one reason why it is widely 
used. However, the physical basis for this model is some- 
what murky. By contrast, fossil electrons do not require any 
new (unknown) acceleration physics; in fact, they mttst exist, 
and the cooling physics which governs their post-injection 
evolution is well established. 

Although it fares much better than other models, our 
fiducial reacceleration model still appears to somewhat un- 
derpredict the luminosities of relics at low Mach number. 
The observations lie well within model uncertainties; besides 
possible variations in area, temperature, density and B-field 
in the downstream plasma, we have not taken into account 
the important effects of variations in relic size along the line 
of sight, viewing geometry, as well as the fact that the fos- 



In reality, of course, much of the change comes from the other 
(uncertain) degenerate parameters. 
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sil electron abundance can vary by a factor of ~ 10 (35T2}- 
Nonetheless, it is somewhat striking that at face value the 
observations do not appear to show any luminosity trend 
with Mach number. Several points are worth noting. Firstly, 
note that many of the relics have spectral indices that are 
spatially unresolved. While we have attempted to correct for 
the effects of cooling, depending on viewing geometry and 
other complications not in our simple model, spectral steep- 
ening due to cooling could be larger. Thus, the shock Mach 
number could be underestimated. Secondly, this could sim- 
ply be a selection effect, given that flux limits typically trans- 
late to luminosity limits of Pi. 4 Si 10^° ergs~^ Hz~^ (except 
for nearby clusters). Indeed, counterexamples to the ob- 
served bright low Mach number relics also exist. For in- 
stance, in A2146 there is a pair of shocks with Mach numbers 
M = 2.1, 1.6 which have been unambiguously detected in 
X-ray observations, but for which no diffuse radio counter- 

! iart has yet been i dentified in deep GMRT observation^^ 
Russeh et ai]|201ll ). Volume-limited samples or shocks se- 
lected in X-ray rather than radio would be required to pin 
down selection bias. Thirdly, characterizing shocks by a sin- 
gle Mach number might be too simplistic; in reality shock 
fronts are curved wi th a spatially varying Mach number 
jSkillman et al.ll2012l ). Finally, we shall see below that the 
non-linear mapping between Ml and Pi. 4 imply that despite 
the predominance of low Mach number shocks, most relics in 
fact cluster about a limiting luminosity close to the asymp- 
totic value of Pi. 4. 

We now make some approximate calculations for the 
radio relic luminosity function. The differential number den- 
sity ofj;eli^sas_a_function of Mach number can be written 
as ||Rvu et al.ll20o3 ): 

i!L = ±^ (37) 
dX *s dM ' ^ ' 

where S is the total shock surface area divided by the volume 
of the simulation box. Thus, 1/S has units of length and can 
be thought of as the mean separation between shocks. We 
Eissume a typical area for relic shocks of 'l/g ~ 2Mpc^. We 
made a fit for the Mach number dependent differential shock 
surface. It is derived from the shock surface of gas with a pre- 
shock temperature of Ti ^ lO^K in a box of size 85 h'^Mpc 
l|Kang fc Rvull201ll '): 

The differential number density of radio relics as a function 
of luminosity is: 

dn dn dj\4 



dPi., 



dM dPi. 



(39) 



where dA^/dPi.4 depends on the acceleration model; we use 
the results of Fig. [15] for reacceleration and direct injection 
('M-aR'). The total number of radio relics inside a cosmo- 
logical box of (1 Gpc)'^ and with a luminosity above Pi. 4 



A^(> P1.4) = (1 Gpcf X J dPi.. 



dn 
dPT^ 



(40) 



The results are shown in Fig. [TB] Two features stand out: (1) 

^■^ Note that the absence of radio emission given the flux limits 
is consistent with our model uncertainties. 
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Figure 15. Luminosity from radio relics as a function of Mach 
number, for reacceleration (solid lines), and direct injection with 
the 'M-aR' (dotted lines) and 'M-const' model (dashed lines). We 
also show observations where the intrinsic spectral index can be 
accurately determined, as in Table [2] (flUed diamonds; labelled), 
and those where we correct for the effects of cooling (open di- 
amonds; upper limits). To bracket the uncertainty in parame- 
ters we vary parameters (lower limit, fiducial value, upper limit) 
and color code the resulting intervals for our reacceleration sce- 
nario of fossils (orange), direct injection scenario (yellow) and 
overlap region (red); shock area 

*S = (0.3, 2,4)Mpc2, down- 
stream temperature kT2 = (1, 5, 10) keV, electron number den- 
sity ric = (0.1, 1, 5) X 10^** cm^^, and B2 = (1, 5, 7) nG (magnetic 
field). 



the reacceleration and direct injection models do not differ 
significantly in shape, even at the faint end. (2) The model 
is strongly discrepant with observations for low luminosities. 
Let us discuss these in turn. 

How can we understand the shape of the relic luminos- 
ity function? Its features can be broadly understood from 
the results of Fig . 1151 C onsider the crosses on the curves for 
the iKang fc Rvul (|201lh parametrization in Fig. 1161 These 
mark Mach numbers in unit intervals, starting from 2 (3) 
for reacceleration (direct injection). The convex shape of 
the luminosity function is due to the fact that luminosity 
falls sharply at low Mach number (even for the reaccelera- 
tion model). Thus, a narrow range in Mach number corre- 
sponds to a wide range in luminosity, which 'stretches out' 
N(> Pi. 4) to the left at low luminosity and leads to its fiat 
slope, since there are relatively few clusters in an interval 
dL. This amount of 'stretching' produced by this non-linear 
transformation differs between the two models, as evidenced 
by the different Mach number ticks. However, the small dif- 
ference in faint end slopes this produces is insufficient to 
serve as a robust discriminant between models. The offset 
on the X-axis, due to different limiting luminosities, simply 
reflects assumptions about asymptotic acceleration efficien- 
cies. The differential luminosity function is similarly poor 
at distinguishing between models, where the main feature 
is a sharp peak at the limiting luminosity. The bottom line 
is that the non-linear mapping between M and Pi. 4 causes 
most luminosities to cluster about a characteristic value, 
which mitigates the efficacy of the luminosity function as a 
discriminant between models. 
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The discrepancy of our model predictions with obser- 
vations is interesting. Such over-prediction of relics also af- 
flicts most current models of the relic luminosity functiorF^. 
Part of the reason may well be that current surveys are 
simply highly incomplete, particularly at low luminosities. 
Note also that: (1) the AS/ AM function diff ers signifi- 
cantly between paramet r izations (compare e.g. Ryu et ahl 
I2OO3I : iKang fc Rvull201ll : lArava-Melo et al.ll2012l ). perhaps 
due to differences in the simulation method, shock finding 
algorithms, and resolution. Given these differences between 
simulations, it is perhaps not surprising that they also dis- 
agree with observations. Also note that our median distribu- 
tion function is valid only for the cluster outskirts (which is 
not singled out by the temperature cut for AS /AM); the dis- 
tribution function closer in suffers increased Coulomb cool- 
ing and has a lower normalization. (2) We have assumed 
fixed values for most relic parameters, whereas in reality 
they should have broad distributions, which would produce 
a tail at high luminosities and smooth out the steep number 
density function. The parameters we adopt are likely biased 
by selection effects. For instance, the magnetic field could be 
highly patchy and inhomogeneous, and only regions which 
have a strong pre-existing field (to be amplified as a shock) 
are visible as radio relics. 

The above issues are well beyond the scope of this 
paper, though they deserve careful study in the future. 
Our main conclusion is that spectrally resolved observa- 
tions which constrain Mach numbers can distinguish be- 
tween reacceleration and direct injection, while the relic lu- 
minosity function, even with outstanding statistics, cannot. 



7 DISCUSSION: INSIGHTS FROM 
HELIOSPHERIC OBSERVATIONS 

Here, we discuss some uncertainties in our adopted shock 
physics and numerical method. In particular, we consider 
our assumptions that the acceleration efficiency plummets at 
weak shocks, and that the shock surface can be characterized 
by a single Mach number with little small-scale spatial or 
temporal structure. We play particular attention to insights 
gained from heliospheric observations. 



7.1 Acceleration Efficiency at Weak Shocks 

Are physical processes in the solar system relevant to the 
ICM? At first blush, the two may appear to differ signifi- 
cantly in their magnetization. While the plasma in the solar 
wind is strongly magnetized with a plasma beta parame- 
ter of /? = Pth/Ps ~ 0.1 — 1 (i.e., t he ratio of thermal- 
to-magnetic pressure, lTreumannll2009l '). the bulk of the in- 



As can also be inferred from the fact that we obtain 
roughly similar numbers. For instance, the number of relics 
with a luminosity Pi. 4 > 10^*^ ergs~^ Hz~^, within a volume 
of (500)'^ (fe~^Mpc)'^, is about 10 00, consistent with what was 
found in lArava-Melo et al] l l2012f) for an acceleration efficiency 
of Ce ~ 0.1 per cent. For LOFAR Tier 1 that has a high sen - 
sitivity of about 0.5 mjy at 240 MHz llMorganti et al.ll201(il) 
(~ 5 X 10^* erg s~^ for an average source distance of 300 Mpc), we 
expect about 3000 radio relics could be visibl e within a volume of 
(500)3(?t-^Mpc)3, compared with ~ 2500 in lNuza et al] l l201lf) . 



observations 
reacc. fossil 
direct inj. 
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Figure 16. Cumulative number density of radio relics in a 
l{h~^ Gpc)'' box. We use the differential shock surface model 
from Kang and Ryu (equation I38I I. with reaccelerated fossil elec- 
trons and direct injection ('M-aR' model) and parameters typical 
for a relic (see Fig. I15I I. Crosses indicate Mach numbers at unit 
intervals, starting with the numbers on the left. For comparison, 
we show the NVSS relic luminosity function (+). 



tra cluster plasma is only weakly magnetized with P > 50 
as infe rred from Faraday ro tation measure studies (e.g., in 
Coma, iBonafede et~alll2010t ). Those, however, are biased to- 
wards the denser core regions and may not be represen- 
tative of cluster outskirts where relics form. In fact, the 
plasma beta parameters at the relics' position are typi- 
cally an order of magnitude smaller, /3 ~ 5. For instance, 
comparing the limit on the inverse Compton X-ray emis- 
sion with the measured radio synchrotron emission for the 
northwest radio relic in A3667 yields a magnetic field es- 
timate of B > S fiG (or equivalen tly P < 3 downstream 
the shock; iFinoguenov et al] |2010| ). Equating the advec- 
tion and synchrotron cooling timescales while deproject- 
ing the relic width for the northern relic in the cluster 
CIZA J224 2.8-F5301 yields -B ~ 5 - 7 fiG (or equivalently 
P ~ 6 ± 3; van Weeren et al] |2010| : lAkamatsu fc Kawaharal 
I2OIII : lOgrean e t al. 2013^^ 

What do heliospheric observations tell us about shock 
acceleration efficiency as a function of Mach number? Ex- 
tensive spacecraft observations in the solar wind have es- 
tablished that electron acceleration is effectively absent at 
quast-parallel weak shocks, i.e. when the upstream magnetic 
field is approximately paral lel to the shock surface nor- 



mal (ISarris fc Krimigid Il985l : ICosling et al] 1 19891: iKrimigis 



199d: lOka et al.ll200d : see a so a recent review bv lTreumannI 
|200£ ) . Electron heating is mainly dominated by conservation 
of magnetic moment (Feldman 1985). Quast-perpendicular 
weak shocks (with Alvenic Mach numbers Ma ^ 3 — 4) 



show electron acceleration as a result of either shock drift 
or surfing accel eration, but diffusive shock accelerat i on ap - 
pears unlikely (| S arris fc KrimigisI [l985l : iTreumannI bOOgj ). 
This limits the maximum electron energy to that achievable 
by shock drift or surfi ng acceleration, wh ich also depends on 
the shock ramp width (|Lever et al.ll200ll ). Thus, weak shocks 
in the solar system do not result in efficient electron accel- 
eration (through DSA), consistent with our assumptions for 
the ICM. One possible interpretation of this is that in order 
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to solve the injection problem of electrons for the DSA mech- 
anism, some process operating under quasi-perpendicular 
conditions (such as shock drift or surfing acceleration) would 
be necessary. By contrast, (rare) strong shocks in the so- 
lar system show significant local electron acceleration up to 
particle energies of a few MeV (see GEOTAIL observations 
of an inte rplanetary shock wit h an Alfenic Mach number 
Ma ~ 6. IShimada et all Il999l '). even under quasi-parallel 
magnetic conditions (as demonstrated by C assini observa- 
tions of Saturn's bow shock with Ma ~ 100, [Masters et al.l 
I2OI3I ). This may hint at diffusive shock acceleration as the 
underlying mechan ism, in agreement with predictions by 
theoretical studies ([Amano fc Hoshindlioiol ). 

These observations suggest that magnetic geometry 
strongly impacts acceleration efficiencies — at least locally 
at the shock where heliospheric observations are providing 
measurements. In Appendix [B] we illustrate how a small 
variation in a geometric parameter can significantly change 
the distribution function of fossil electrons. Nonetheless, we 
do not propagate uncertainties in magnetic geometry across 
our calculations, and only consider DSA as a source of accel- 
eration. We believe this is justifiable in observed radio relics, 
on two grounds: (1) The radio spectral slopes are consistent 
with power-law electron spectra due to diffusive shock ac- 
celeration. The high Lorentz factors of 7 ~ 10* required by 
the radio data cannot be achieved by shock drift or surfing 
acceleration. (2) The smoothness of the radio intensity and 
polarization along the shock surface on scales of ~ 2 Mpc 
argues for a robust acceleration process, independent of the 
local field geometry. Nonetheless, magnetic geometry is an 
important caveat in considering the ensemble of possible ra- 
dio relics: for instance, perhaps only a fraction with favorable 
geometry are visible. 

7.2 Shock Surface 

In this paper, we have characterized shocks by a single Mach 
number (this is particularly important in [JB] when we com- 
pare our model to observations) . Ev en in the DSA frame- 
work, this is an over-simplification (jSkillman et al.ll2012l ). 
More generally, the shock surface is likely to be complex 
and dynamic. For instance, the crossing of the termination 
shock in the solar wind by Voyager 2 revealed the presence 
of a highly dynamic shock surface. Instead of a single cross- 
ing of a stationary shock, Voyager 2 passed the shock sev- 
eral times on its ballistic orbit into the heliosheath. The 
data reveal a complex, rippled, quasi-perpendicular super- 
critical magnetohydrodynamic shock that undergoe s refor- 
mation on a scale of a few hours (|Burlaga et al ]|2008l l. which 
corresponds to a characteristic reformation length scale of 
~ 10^^ cm ~ 0.01 AU (for a solar wind shock velocity of 
300 km s~^). Supercritical shocks are not capable of gener- 
ating sufficient dissipation for the retardation, thermaliza- 
tion and entropy increase associated with a shock transi- 
tion. When the flow exceeds a critical Mach number for a 
given magnetic obliquity, the supercritical shock starts re- 
flecting particles back upstream. Beyond a second critical 
Mach number, whistlers accumulate at the shock front and 

A 'supercritical' shock is one in which gyrating ions represent 
the primary dissipation process. 



reform periodically (|Treumannll2009l ). In addition to these 
microscopic shock reformation mechanisms, there should be 
macroscopic reformation because of the varying ram pres- 
sure and magnetic field fluctuations in the upstream that 
cause it to reform on gradient length scales of the infalling 
structures, which likely exhibits a characteristic size spec- 
trum. These are complications well beyond the scope of this 
paper, but they illustrate possible complexities associated 
with the shock surface. 

More mundanely, our results can depend on the numer- 
ical method used to track the injection of CRs. There are 
known differences in the morphology of accretion shocks in 
SPH versus mesh-based finite volume methods. Entropy pro- 
files show a sharp peak located around ~ 2 — 3 i?vir in mesh- 
based techniques, while the profiles in SPH simulations — 
although similar in shape — are smoothed across a sizable 
larger volume, and the volume-weighted Mach number dis- 
tribution, particularly for external shock s in low density re - 
gions, have different trends in each code l|Vazza et al.ll20l"ll V 
However, the Mach number distribution weighted by dissi- 
pated energy (which is the most important quantity here) 
is reassuringly similar for those different numerical schemes 
owing to the conservativ e nature of the implemented equa- 
tions l|Vazza et al.ll201ll ). For the resolution used here, we 
hence do not expect dramatic changes, though direct confir- 
mation is of course preferable. 



8 CONCLUSIONS 

In this paper, we consider giant radio relics in galaxy clus- 
ters as the manifestation of reaccelerated fossil relativistic 
electrons, rather than direct injection from the thermal pool. 
This idea is not new a nd has already been proposed in the lit- 
eraturc several times (Ensslin et al."l998': 'Markevitc h et al. 
2005; Giacintucci ct al. 2008; Kang & Rvu 2011; Kan g et al. 
2012l :l van Weeren et al.ll2012l . e.g.,). What has been conspic- 
uously missing from the literature is an actual calculation 
of the fossil electron distribution function, and whether it 
has the correct shape and amplitude to account for the ob- 
servations. This paper represents the first attempt to do so. 
We run a suite of cosmological simulations where cosmic 
ray electron injection and radiative and Coulomb cooling 
are tracked in post-processing. 

Our principal findings are as follows: 

• The fossil CRe population is substantial. Without any 
fine-tuning, it is at the right level to explain observations. 

• Reaccelerated fossil electrons are competitive with di- 
rect injection at high Mach numbers, and overwhelmingly 
dominate by several orders of magnitude at low Mach num- 
bers M ^ 3. Because it is a relative comparison, this con- 
clusion is fairly robust to the (potentially uncertain) accel- 
eration model, and depends on well-understood structure 
formation and cooling physics. Given that low Mach num- 
ber shocks are strongly prevalent in clusters, we predict that 
LOFAR should find many more unexpectedly bright steep- 
spectrum radio relics, which cannot be explained by direct 
injection. 

• The fossil distribution function has a generic shape 
with distinct regimes, where different cooling processes dom- 
inate: (1) sub-relativistic Coulomb cooling (p ^ 1), (2) rela- 
tivistic Coulomb cooling (l^p^lO), (3) quasi-adiabatic 
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(lO^ip^lOO), (4) inverse Compton/synchrotron cooling 
(p ^ 100). Although fossil electrons peak in number den- 
sity at p ~ 1, due to the sharp change in the nature of 
Coulomb cooling, electrons at higher energies (1 ^ p ^ 100 
and 10 ^ p ^ 10* for strong and weak shocks respectively) 
contribute most to the observationally relevant portion of 
the reaccelerated spectrum. Cooling times in these regimes 
are relatively long, and thus it is possible to use simulation 
outputs with relatively poor time resolution (At = 100 Myr, 
1 Gyr). 

• Since both injection and Coulomb cooling change on the 
(relatively long) dynamical timescale, they can be approxi- 
mated as steady. We compare a self-similar analytic model 
where this approximation is made and find excellent agree- 
ment with the simulations. This enables extremely rapid es- 
timates. Regime (1), (2) & (4) are in steady-state, while (3) 
grows monotonically with time. 

• There can be up to a factor of ~ 10 scatter in fossil 
electron abundance, depending on a cluster's merger history; 
cool core clusters have a lower abundance. Spatial variations 
within a cluster can also be significant. 

Given these considerations, we strongly advocate that 
fossil electrons be considered a key ingredient in interpreting 
radio relic observations, rather than an exotic afterthought. 
Failure to account for the fossil electrons will lead to er- 
roneous conclusions about the nature of particle accelera- 
tion at weak shocks. Several extensions immediately sug- 
gest themselves. Seed electrons are needed for turbulent 
reacceleration, as in a pr ominent model for radio halos 
(|Brunetti fc Lazarianll2007l ): we are presently studying if fos- 
sils from structure formation are adequate. Maps, more de- 
tailed and realistic luminosity functions and observational 
predictions are needed. As discussed in H7.21 confirmation 
and extension of these results with a grid code would be 
most welcome. Our main hope is to stimulate others to pur- 
sue such lines of inquiry. 
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APPENDIX A: COOLING 

In this section we explain our implementation of cooling 
in the "impulsive" scenario in more detail. These are used 
to evolve the CRe population in our simulations with finite 
time resolution. We develop analytic expressions which allow 
very rapid calculations. 



Al Coulomb cooling 

We start by deriving the shift in momentum pi at time ti to 
momentum pf at time tf due to Coulomb cooling by inte- 
grating equation ((23]): 



dp' 
fec(p') 



dtbc{t) = F{tf,U) , 



(Al) 



where bcit) — ne{t)/n, and n an auxiliary variable intro- 
duced to have the correct units. The momentum dependent 
part is determined by 



_dp_ 

bc(p') 



nflx dp' 
bc{p',t) 



(A2) 



* cm and 



for a fixed electron number density wax = 10 
where 6c (p, t) is given by equation (I24|l . The integral over 
time is performed over snapshots that are discrete in time. 
Hence, we approximate the time integral for Coulomb cool- 
ing between time ti and by a discrete sum: 



F{tf,U)^ J2 At,[focfe-i) + 6cfe)]/2, 



(A3) 



where j denotes the summation index that run over all snap- 
shots between ti and tf. Here we have approximated the time 
integrated cooling rate between two snapshots by the mean. 

First we look at the momentum integral, which we ap- 
proximate with a simple fit accurate to 20 per cent within 
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10"^ <p< lO"*, and 10"'^ < < 10' 



/ 



dp' 

bc{p') 



where , 



7.73 X 10" s, and (ji = i.Q . 



(A4) 



We introduce a characteristic momentum for the Coulomb 
cooling, pcoui, such that 



PCoul 



= F{tf,U) , where (A5) 




and 



K = <pF^a.nd J = 27 (ji + A . 



(A7) 



Note that pcoui ~ bc{p,t)At in the relativistic regime. The 
shift in momentum is then given by 



Apcoui 
Pi 



Pf -pi 



(A8) 



6 

+ 2F' 



22/3 {j'p'_2 F'^ + f'^^K' J'^ 
2^F' 



1/3 



1 + 



J' F' -2F''-^ + 2,2y/TK' 



and 

F' = 
K' = 



PCoul 



(tl + pj (l> + p^^^l 

(AF'^and j' = 27 <?!) + 4 



(A9) 

(AlO) 
(All) 



A2 Inverse Compton Cooling 

We similarly derive the shift in momentum due to inverse 
Compton cooling by integrating equation p6|l from the red- 
shift Zi = z{ti), where the electrons are injected, to a later 
time Zf = z(tf). We define a characteristic momentum pic, 
where: 



Pic ~ Pf Pi^ Hopj (1 



{l + z)*dz 



zj (1 + z) ^^m{1 + zf + VLf, 

(A12) 

Note that pic ~ l/(6ic,o At) for relatively short timescales 
(where At < 1 Gyr, and tiic,o is given by equation (|25p '). 
After the time At = tf — ti has elapsed, all CRes with a 
momentum above pic have cooled to lower momentum. The 

_ 2 

shift in momentum is given by Apic ~ — ^ — . 



A3 Cooled Distribution Function 

Given an initial energy pi of an electron at time ti, equa- 
tion (|23|) can be integrated to give the value of pf at a later 
time tf. The differential population density for relativistic 




/',c(t) 



Figure Al. Spectral behavior of the cooled CR electron spec- 
tra for different momentum regions. There are several distinct 
momentum regimes: (1) the low momentum regime (p < 1) is 
dominated by sub-relativistic Coulomb losses, which cool the CR 
electrons very efficiently, (2) the less efficient relativistic Coulomb 
cooling regime (1 ^ P ^ 10), (3) the adiabatic regime (10 < 
P ^ 10^), which preserves the injected distribution function, and 
where Coulomb and inverse Compton/synchrotron losses are less 
important. (4) the inverse Compton/synchrotron cooling regime 
(p > lO'^). We assume a power-law index of the injected CR elec- 
trons of a = 2.5. 



electrons is then given by 

/inj,CRe(Pf,t/,tO = /inj,CRc(Pi,tO 



dpf 



(A13) 



where 

/inj,CRc(Pi, ti) = /inj,CH.o(pf - ApiC - ApCoul , U) , (A14) 

and 

dpi 



dpi 



P\c 



14- 



j (pic - Pf) 

^ 3<7i-fp| (<?i+Pi(pcoui)')' 



Pf 



Pi(pcoui)/ ((^ + pf)^ 3(?!) +pi(pcoui)^ 



(A15) 



The total electron spectrum is derived from the sum of all 
individually cooled injected spectra, starting from the time 
of injection ti until a later time tf, 



/CRe(pf , i/) = X] /inj,CRe(pf , t/, tj 



(A16) 



APPENDIX B: 
PARAMETER 



DEPENDENCE ON INJECTION 



Here we investigate the importance of the parameter that 
controls the leakage of the injection process, ea = Bo /B± ■ 
It is given by the ratio of the amplitude of the downstream 
MHD wave turbulence, Bo, to the magnetic field along the 
shock normal, B^. The physical range of this parameter 
is quite uncertain due to complex plasma interactions. We 
adopt a value of eb = 0.23 in the paper, which corresponds 
to about a maximum acceleration efficiency of electrons of 
Ce,inj = 0.01. In Fig. IBll we show how a small variation in 
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Figure Bl. The CR electron dependence on the acceleration 
parameter eb = / . We show the median fossil CR electron 
spectra with full cooling for a large post-merging cluster at z = 
in the region between (1.3-1.5)-R200 ■ The spectra are shown for 
different values of cb; the value used in the paper eb = 0.23 
(green dotted line), ea = 0.22 (blue solid line), and eb = 0.24 
(red dashed line). 
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this parameter change the distribution function of the fos- 
sil CRes. We find a strong dependence on this parameter. 
However, since eb is essentially degenerate with Cc,mj, it can- 
not vary too much from our assumed value, as that would 
confiict with observations of supernova remnants and radio 
relics. 



APPENDIX C: ANALYTIC FIT TO CR 
INJECTION RATE 

In this Appendix we present an analytic formula for the 
time-averaged (over the entire cluster history) CRe injection 
rate. In the left panel of Fig. ICll we show the injection rate 
of each cluster as a function of momentum, where we have 
normalized all curves to have the same value at p = 10 
(adiabatic regime). There is remarkably little difference in 
the spectral shape, indicating that all clusters have a similar 
mix of strong and weak shocks. We use a triple power-law 
fit to capture the spectral shape of the median injection rate 
over the cluster sample: 

(3mj,ad(P,t/,i»)/(3mj,ad(10,t/,ti) = "^Aip'"' , 

i 

where 

Ao = 194.4 , Ai = 172.8 , A2 = 68.4 , and 
ao = 2.35, ai = 2.5, aa = 2.7 . (CI) 

This fits to better than 5 per cent over the entire depicted 
momentum range. 

The injection rate at momentum p = 10 where 
the electron number density is renormalized, Qmj,ad(p = 
10, t/, ii)(nc,norm/".c, cluster), is showu in the right panel of 
Fig. ICll There is no clear trend with cluster mass or mor- 
phology. Hence, we adopt a median value of (3inj,ad(p = 
10, t/,ii)(?ic,norm/nc,ciustor(t/)) = 5 X 10~^° cm"^ s~\ whcrc 
we renormalized the density with ne,norm = 10~^ cm""^. 



© 2011 RAS, MNRAS 000.[Tlf23l 



24 A. Pinzke, S.P. Oh and C. Pfrommer 



10' 

0.10 F 

0.05 h 
0.00 i- 
-0.05 h 

-o.ioL 
10-- 



M-aR 



h simulation - median : 
simulation - individual - 
fit 

= 2.70 1 
= 2.50 = 
= 2.35 




10--' 



10' 



10' 
: P, I m c 



10' 



10* 



V 












simulation: 


■ M-aR 








+ cool core clusters - 










X merging clusters 










X 






X X 


X 


+ 

1 

++ 






+ X 








X 




X 







10" 



10-' 10" 
[10" Mel 



Figure CI. Cosmic ray electron injection rate. We show the median injection rate of CR electrons in the cluster region between 
(1.3-1.5)i?200 for our sample of 14 simulated galaxy clusters. The upper panel in the left figure shows the spectra as a function 
of dimensionless CR electron momenta. The injection rate for each cluster is normalized at p = 1 (grey solid lines). The black 
crosses show the median of the normalized injection rate across our cluster sample. The blue line shows the best fit triple power- 
law to the median spectrum, where the red, orange and green lines show each of the fitted power-laws. The bottom left panel 
shows the difference between the relative fit and the simulation data (blue solid line) which amounts to less than 5 percent. The 



right figure shows the injection rate, < Qinj,ad(p = Iji/i^i) >timG (f^e 



e, norm /'^e, cluster 



(tf)), renormalized with the electron density 



'^c,iiorm = 10~ cm~'^, and averaged over time, as a function of galaxy cluster mass M200 a^nd cluster morphology. The injection 
rate for each merging cluster (red crosses) and each cool core cluster (blue -|-) are shown. Note that our limited statistical sample 
of clusters shows no clear trend with mass. 
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